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Abstract. The Bogomolov Conjecture is a finiteness statement about algebraic points of 
small height on a smooth complete curve defined over a global field. We verify an effective 
form of the Bogomolov Conjecture for all curves of genus at most 4 over a function field 
of characteristic zero. We recover the known result for genus 2 curves and in many cases 
improve upon the known bound for genus 3 curves. For many curves of genus 4 with bad 
reduction, the conjecture was previously unproved. 



1. Introduction 

1.1. The Conjecture and Main Theorem. Fix an algebraically closed field k of charac- 
teristic zero and a smooth proper connected curve Y/k. Define K to be the field of rational 
functions on Y. Let C be a smooth proper geometrically connected curve of genus at least 
2 over the function field K . Choose a divisor D of degree 1 on C = C x K K and con- 
sider the embedding of C into its Jacobian Jac(C) = Pic°(C) given on geometric points by 
Jd{x) = [x] — D. Define 

a'(D) = liminf h (Jd{x)) , 
xec(K) 

where h is the canonical Neron-Tate height on the Jacobian associated to the symmetric 
ample divisor + [— 1]*@. As C(K) may not be countable, the liminf is taken to mean the 
limit over the directed set of all cofinite subsets of C(K) of the infimum of the heights of 
points in such a subset. Recall that C is called constant if there is a curve Co defined over 
the constant field k and a finite extension K'/K such that Ck 1 = Co x /c K' . Write Div 1 (C) 
for the set of divisors of degree 1 on C. We wish to investigate the 

Geometric Bogomolov Conjecture ([3]). If C is not a constant curve, then 

inf _ a'(D) > 0. 

In response, we have proved the following result: 

Theorem 1.1. Let K be a function field of characteristic zero. Then the Geometric Bo- 
gomolov Conjecture is true for all curves C/K of genus 2 < g < 4. Moreover, if C is not 
a constant curve, there is an effectively computable positive lower bound for a'(D) that is 
uniform in D. 



A more precise form of this result is given by Theorem 1.2 below. 
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It is worth noting that the dependence on D in the Geometric Bogomolov Conjecture is 
superficial. Indeed, it follows from [17, Thm. 5.6] that for any degree-1 divisor D, 

a\D)>-a' ( Kr 



2 \2g-2 

where Kq is a canonical divisor on C and Ko/(2g — 2) is a degree-1 divisor with rational 
coefficients. Thus is suffices to obtain a positive lower bound for the single divisor £ = 
Kc/(2g — 2). Now we can reformulate the Geometric Bogomolov Conjecture in the following 
more intuitive form: 

Height Gap Principle. If C is not a constant curve and £ = Kc/(2g — 2), then there 
exists e > such that for any x G C(K), 

htit(x))^0=>htit(x))>e. 



The Height Gap Principle is true for all curves of genus at most 4 by Theorem 1.1 
However, the theorem does not provide any information on the algebraic point of smallest 
positive height, so we cannot effectively choose e in the statement of the Height Gap Principle. 

1.2. A More Precise Statement and Some History. We now provide a more precise 



statement of Theorem IT We continue to use the notation from the previous section, but let 
us stress that we still assume the curve C has genus at least 2. By the Semistable Reduction 
Theorem we may pass to a finite extension field K' over which Ck> = C x K K' has semistable 
reduction. One can do this effectively by choosing K' so that all of the 12-torsion points of 
Jac(C) are rational. For references on semistable reduction theory, including proofs of these 
facts, see [13], expose n° 1] and [8, §9.3.3, 10.3, 10.4]. 

Let Y'/k be a smooth proper curve with field of rational functions K' . To say that Ck> has 
semistable reduction means there is a projective surface X'/k and a proper flat morphism 
/ : X' -> Y' so that 

• / has generic fiber isomorphic to Ck<- 

• The fibers of / are connected and reduced with only nodal singularities. 

• If Z is an irreducible component of a fiber and Z = P 1 , then Z meets the other 
components of the fiber in at least 2 points. 

If we assume further that X' is a smooth surface (over k), then such a morphism / : X' — > Y' 
is unique up to canonical isomorphism, and it may be characterized as the minimal regular 
model of Ck> over Y'. 

We divide the fiber singularities of / into different types as follows. Choose a point 
y G Y'(k). The partial normalization of the fiber f~ x {y) at a node p is the fc-scheme 
given by resolving the singularity at p. We say p is of type if the partial normalization 
at p is connected. Otherwise, the partial normalization at p has two connected components. 
If the fiber f~ 1 (y) has arithmetic genus g, then one component of the partial normalization 
has arithmetic genus % and the other has arithmetic genus g — i. We may assume i < g — i. 
In this case we say the node p is of type i. Let 5i = 5i(X'/Y') be the total number of 
nodes of type i in all fibers. By uniqueness of the minimal regular model, the numbers Si 
are well-defined invariants of Ck>- 
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Theorem 1.2. Let K be a function field of characteristic zero over the algebraically closed 
constant field k. Let C/K be a smooth proper geometrically connected curve of genus g 
with 2 < g < 4. Suppose that K' is a finite extension of K over which Ck> has semistable 
reduction. Define Y' to be the smooth curve with function field K' , f : X' —* Y' the minimal 
regular (semistable) model of Ck>, and d = [K' : K\. Then 

3 

d(g - 

1 I . ^ 2i{g-i) 



if f is smooth 



inf _ a\D) > < 

DeDiv^c) 



2d(2g + I) 



c(g)5o + 5i J unconditionally 

i6(0, ff /2] 



with c(2) = c(3) = and c(4) = ^. In particular, the Geometric Bogomolov Conjecture 
is true for curves of genus at most 4. 

Results of this type have been known for about a dozen years now, and the novelty of the 
present paper is twofold. First, it gives an algorithm for verifying the Geometric Bogomolov 
Conjecture for all curves of a fixed genus. Second, the method given here seems to admit 
a generalization, although we do not yet completely understand it. We are in a position 
to make the following strong effectivity conjecture, which we do not claim is in any way 
optimal: 

Conjecture 1.3. Let C/K be a smooth proper geometrically connected curve of genus g > 2 
that admits semistable reduction over K . Then 

inf a'(D)> , 1 J 9 —6 0+ V 2% - %) 5 % 
o^ic) J -2(2, + l) ^27, ^ g 

Now let us summarize the previous results relevant to the conjecture. For simplicity, let 
us assume that C has semistable reduction over K (so that d = 1 in the theorem) and that 
/ : X — > Y is the minimal regular model. 

(a) (Parsin, [12J) For any g > 2, if / is smooth then the relative dualizing sheaf ujx/y is 
ample and its self-intersection number satisfies ojx/y-^x/y > 12. By [T71 Thm. 5.6] 
this implies 

inf _ a\D) > 



(See also Theorem 3.2 below and the subsequent discussion.) 



(b) (Moriwaki, [9]) If the genus of C is 2, then / is not smooth and 

inf a'(D) > -L (J + Sl ) . 

(c) (Moriwaki, [TU1 [TT]) If the dual graph of each closed fiber of the stable model of / 
consists of a tree with loop edges attached, then 

[9/2] 



inf a\D)> , 1 A g —5o + Y 2 -^> 



neDiv^c) 2(2flr + l) I 6g g 




(d) (Yamaki, [16]) If C is hyperelliptic of genus g > 3 and the hyperelliptic involution i 
extends to the family /, then 



inf _ a'{D) > ' 



DeDiv 1 ^) 4g(2g + l) 



2g-5 



[(fl-i)/2] 



£o + £ {(<7-i-i)-i}& 



[9/2] 

- a(5){^ - Co} + 4i ( 5 ~ 0* 
i=i 

where £ = £o{X/Y) is the number of nodes of type that are fixed by the involution 
and £j = £j(X/Y) is the number of pairs of nodes {p, t(p)} of type not fixed by the 
involution and such that the partial normalization at {p, i(p)} yields a curve of genus 
j and a curve of genus g — 1 — j. The constant a(g) is to be interpreted as when 
g = 3, 4 and as ^p- when g > 5. (This follows from Yamaki's formulas upon noting 

that 5 - £ = Ei=T 1)/2] 6i, by definition.) 

(e) (Yamaki, [T5]) If the genus of C is 3 and C is not hyperelliptic, then 

inf a '(D)>-U + k. 
DeDiv^c) 792 6 

(f) (Gubler, [7]) If the Jacobian of C has totally degenerate reduction over y e Y(k), 
then the Geometric Bogomolov Conjecture holds for C. This hypothesis is equivalent 
to saying the first Betti number of the reduction graph equals g, the genus of Cf] 

For genus 2 curves, we obtain the same coefficient on 5q as Moriwaki. This can be explained 
by the fact that every curve of genus 2 is hyperelliptic, and one of the key inequalities used 
to prove our result becomes an equality for hyperelliptic curves. (It is the second half of 



inequality fll-iP -) For genus 3 curves, we obtain a lower bound of -^5q + which is in 



general neither stronger nor weaker than Yamaki's result. One can recover the result (c) 
above using the computation in [181 Prop. 4.4.3] and the method of this paper. 

Remark 1.4. In all but one of the related articles of Moriwaki and both of the articles of 
Yamaki, they use a slightly different embedding of the curve into its Jacobian and a different 
measure of finiteness than we do. At first glance, the results in their papers will look quite 
different than as stated here. 



1.3. Idea of the Proof. Now we sketch the proof of Theorem |1.2| in order to motivate the 
layout of the paper. The detailed proof will be given in §|3} 

Assume for simplicity that C has semistable reduction over K, and let / : X —>■ Y be 
the minimal regular model of C. If / : X — > Y is smooth, then the result is essentially due 

1 This equivalence does not seem to be written anywhere in the literature, so we give a brief sketch of its 
proof at the request of the referee. The genus of C agrees with the genus of its reduction graph at y, which 
is the first Betti number of the graph plus the sum of the geometric genera of the irreducible components of 
the fiber of X over y. (See ^2.1 and ^3] especially Proposition 3.1 ) To show that the first Betti number of 



the reduction graph equals the genus of C is then equivalent to asserting all of the irreducible components of 
the fiber of X over y are rational. But to say that the Jacobian of C has totally degenerate reduction over 
y is equivalent to saying the abelian rank of the special fiber of the identity component of its Neron model 
is zero, which in turn is equivalent to saying that every irreducible component of the fiber of X over y is 
rational. These facts about the Jacobian can be deduced from [13l expose n° 1, p. 28]. 
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to Parsin, so we may assume that / has at least one singular fiber. We use the following 
inequalities proved by Zhang: 



a '{D) > > , x T V (r„) , (1.1) 

V ' - 4(o - 1) - 2(2o + 1) ^ V W V ; 

where u a is the canonical sheaf of C equipped with the admissible adelic metric in the sense 
of Zhang, {uj a ,u) a ) is the admissible pairing of u> a , and ip (T y ) is a graph invariant associated 
to the reduction graph T y of the fiber f~ 1 {y)- Recall that the reduction graph has vertices 
in bijection with the irreducible components of f~ x (y) and a segment of length 1 between 
two vertices for every node shared by the two corresponding irreducible components. If an 
irreducible component of f~ l {y) has nodal self-intersections, then the corresponding vertex 
of the reduction graph is joined to itself by loops of length 1. For example, T y is a point 
if and only if / is smooth over y, and in this case tp (T y ) =0. As / is generically smooth, 
the above sum over points of Y is finite. To prove Theorem |1.2[ it suffices to produce lower 
bounds for (p (T y ) for all possible (polarized) metric graphs that can arise from fibers of /. 
We state such a lower bound as Theorem 13.41 below. 

The majority of this paper is devoted to exhibiting lower bounds for the invariants <p (T y ) . 



In £2.1-2.3 we gather the definitions and basic facts on metric graphs used to state and prove 
such bounds. In §|3] we give a more complete treatment of the argument sketched above in 
order to pass from a statement about algebraic curves to a statement about graph invariants. 
In ^4 we provide explicit formulas for computing ip (T) that allow one to implement its 
calculation in a computer algebra package. In §5.1| - ^5T3"1 we reduce the proof of the theorem 
to a handful of explicit calculations. These calculations are immensely complicated, albeit 
not very difficult. We perform them in Mathematica [H] and summarize the conclusions 



in §5.4[ Finally, we reproduce the computer code and Mathematica notebooks used for the 
computations. 
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2. Polarized Metric Graphs and Their Invariants 

2.1. Polarized Metric Graphs. In order to fix terminology, we give a very basic intro- 
duction to polarized metric graphs. For a more complete discussion of this topic, see [TJ [2] . 
A metric graph T is a compact connected metric space for which each point p admits a 
neighborhood isometric to one of the form 

U r)V = {te 2nik/v eC:0<t<r,k = l,...,v}, 

for some positive real number r and positive integer v, where we endow U rjV with the path 
metric. The number v = v(p) is called the valence of the point p. For convenience, we also 
allow the topological space consisting of a single point to be called a metric graph; in this 
case, v(p) = 0. 
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When working with metric graphs, it is often easier to "discretize" and work with classical 
combinatorial objects. The following two definitions are slightly non-standard, but they will 
simplify our discussion a great deal. A combinatorial graph is a connected multigraph; 
i.e., a connected graph, possibly with multiple or loop edges. If G is a combinatorial graph, 
we write V(G) and E(G) for the vertex and edge sets, respectively. A weighted graph is a 
combinatorial graph along with a function I : E(G) — > R>o- We call £(e) the length of the 
edge e. 

Any weighted graph G gives rise to a metric graph T first by building a 1-dimensional CW- 
complex using the vertex and edge sets as gluing data. An edge e of G induces a segment 
e C T, which we equip with the Lebesque measure of total mass £(e). The metric on T is 
defined to be the induced path metric. Conversely, given a metric graph T, one obtains a 
weighted graph G as follows. Let V(G) be any nonempty finite subset of T containing all 
points of valence different from 2. Two vertices p\ and p 2 of G are joined by an edge if there 
is a path in T from pi to p 2 that does not pass through any other vertices. Any weighted 
graph constructed in this way from T is called a model of T. 

For any model G of a metric graph T, we have the equality 



This formula can be proved by counting the edges in a spanning tree for G. The number 
b\ = dimjR i? 1 (r, R) is called the first Betti number of T. 

A metric graph is irreducible if it cannot be disconnected (as a topological space) by 
deleting any single point. A metric graph T is cubic if the following two conditions hold: 

• There exists p E T with valence 3. 

• Every p E T has valence equal to 2 or 3. 

This definition implies that a cubic metric graph T admits a 3-regular model in the sense of 
combinatorial graph theory. Indeed, take the vertex set {p E T : v (p) = 3}. (Recall that a 
combinatorial graph is 3-regular if every vertex has valence 3.) 

Let T be a metric graph and choose a model G. Each segment e is isometric to either a 
closed interval or a circle of length £(e). If e is isometric to a segment, a choice of isometry 
e ^ [0,£(e)] gives an orientation on e. We write e~ and e + for the vertices corresponding to 
and £(e), respectively, under the chosen isometry. If e is instead isometric to a circle, we 
define e~ = e + to be the unique marked vertex on e. We define the total length of V to be 
£{Y) = ^2 e £(e). Evidently it does not depend on the choice of model. 

A polarized metric graph, or pm-graph, consists of the data of a metric graph T and 
a vertex weight function q : T — > Z with the following properties: 

• q has finite support in T; 

• q(x) > for all x ET; 

• The canonical divisor K := Xlper (2<?G°) + v (p) — 2) - [p] is effective. 

By way of notation, we will write T = (T, q) for a pm-graph. If Y = {p} is a single point, 
then K is effective if and only if q(p) > 1. Otherwise, it suffices to focus only on points of 
valence 1 to check that K is effective. 
The genus of T is defined to be 



dmiRtf^R) = #E(G) - #V(G) + 1. 



(2.1) 




(2.2) 



per 
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We will abbreviate by g = g (r) when there is no possibility of confusion. It is easy to see 
that the degree of the canonical divisor is deg(K) = 2g — 2. 

Let T = (r, q) be a pm-graph. A model of T is a model of the metric graph T that 
contains the support of q in its vertex set. We will call T irreducible (resp. cubic) if its 
underlying metric graph has this property. 

With T as in the last paragraph, choose a point p G T of valence 2 outside the support 
of the function q. All but finitely many points of T have these two properties. We wish 
to define the type of the point p to be an integer in the interval \0,g /2]. The type of p 



will determine how it contributes to the invariant ip (T). (See Lemma 5.11 ) We say p is 
of type if r \ {p} is connected. Otherwise we may write T = Ti U T^Tor some metric 
subgraphs Ti and T 2 such that Ti fl T 2 = {p}. While the data (Tj,q\r,) does not define a 
pm-graph, the definition of genus ( |2.2[ ) still makes sense. Moreover, one sees that 

9 = 9 (r) = g (Ti, g| r J + g (T 2 , q\r 2 ) , 

so that one of these two integers lies between 1 and g/2, inclusive. Define p to be a point of 
type i, where i is the minimum of g (Ti, glrj and g (T2, q\r 2 )- See Figure [TJfor an example. 



• O • 



Figure 1. Here we have a pm-graph T = (T, q) of genus 4. The indices 1 and 
2 are the values of the function q at those vertices, and we suppose q is zero at 
all other points. A point p in the interior of segment e has type as V \ {p} 
is connected. If p lies in the interior of e' (resp. e"), then it has type 1 (resp. 
type 2) since it splits the graph into a piece of genus 1 and a piece of genus 3 
(resp. two pieces of genus 2). 

Finally, define (T) to be the total length of all points of type i (in the sense of Lebesgue 
measure). Evidently £{T) = Yli^i (X)- (Note that £j (T) for i > is an invariant of the 
pm-graph T, while £0 (T) and £(T) are metric invariants.) 

2.2. Admissible Measures. For this section, fix a pm-graph T = (T,q). We now give a 
very brief description of the theory of admissible measures; for more complete references see 

There is a measure-valued Laplace operator on T. Let / be a continuous function that 
is C 2 outside of some vertex set for T, and such that /" is square-integrable against the 
Lebesgue measure on V. The Laplacian of / is defined as 

A x (f) — —f"(x)dx — ^2a p (f)5 p (x), 

per 

where dx is the Lebesgue measure on T, 5 P is the point mass at p, and <J p (f) is the sum over 
all tangent directions at p of the outward pointing derivatives of /. One of the key properties 
we will use is that A x (f) = if and only if / is a constant function. 

For any signed Borel measure v on T of total mass 1, the Green's function with respect 
to v is the unique function g v : T x T — * K satisfying the following three properties: 
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(a) g v {x,y) is continuous on V x V and symmetric in x and y. 

(b) A x g u (x, y) = 5 y — v for each y e T. 

(c) J r y)dv(x) = for each ?/ G T. 

Here and in what follows it will be convenient to write f{D) = '^2 per n pf(p) if / : T — ► R 
is a function and D = J2 P er n p[p\ * s a divisor on T. As is customary for divisors, we suppose 
n p = for all but finitely many p G T. 

There exist measures whose Green's functions satisfy extra properties. For any divisor 
D on T with deg(D) ^ —2, there exists a unique Borel measure po of total mass 1 and a 
unique constant a = a (T, D\ such that 

9fi D ( x , x ) +9vn( x i D ) = a i xeT. (2.3) 

The measure \xd is called the admissible measure with respect to the divisor D. Of 
particular interest to us is the case D = K, where K is the canonical divisor on V . Note 
that deg K ^ —2 as K is effective for any pm-graph. We will write p = px for the duration 
of this article. 

The effective resistance function r(x,y) — familiar from circuit theory — allows us 
to give an explicit formula for the measure pn- By definition, r(y,z) = j z (y,y), where 
f(x) = j z (x,y) is the fundamental solution to the Laplace equation 

A x f = 5 y - 5 Z , f(z) = 0. 

Let D be a divisor of degree different from —2 and suppose G is a model of T containing the 
support of the divisor D. Let r(e) be the effective resistance between the endpoints of the 
segment e, and let F(e) = 1 — r(e)/£(e). Then we have 

taV ; ^ peV(G) eeE(G) 

In particular, if g = g (r) > 1 is the genus of T and G is a model of T, then 




i{ E ^ +I>(e) gA 

y ^ per ee£(G) v y J 



2.3. The Invariant y?(r). Define three invariants associated to a pm-graph T of genus 
<?>L 

e(r) = JJ r(x,y)5 K (x)dp{y) 
o(r) = ^ ^ r(x,y)dp(x)dp(y) 
<p(F) = 3ga(T)-^{e(T)+£(T)). 

Remark 2.1. We have chosen to use the notation a (T) instead of r (T) as is favored in [18]. 
This is to avoid confusion with the closely related "Tau-invariant" that appears in [21 E]- 



One can show that a (r) agrees with the constant a = a(T,K) appearing in (2.3). Also, 
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the invariant e (r) is the same as the invariant e y appearing in [9] and the invariant e(G, D) 
appearing in [HI [TH1 US], although it requires a small computation to see it. 

Remark 2.2. An important - - albeit trivial — example is the case V = (T,q) in which 
r = {p} consists of a single point. Then T is a pm-graph if and only if q(p) > 1. It follows 
easily that all three of the associated invariants e (T) , a (r) , and (p (r) are zero. 

We now show that a (T) may be removed from the definition of ip (T) . 

Proposition 2.3. For any pm-graph T of genus g, we have 

"I 1 ) = A„<„_^ e I 1 ) 



4007-1) v ' 8g(g-lY 

Proof. Let go : T x T — > W and /io be the admissible Green's function and the admissible 
measure associated to the divisor D = 00 Recall that this means 

(a) g (x,y) is continuous on T x T and symmetric in x and y. 

(b) A x g (x, y) = 8y- /i for each y G V. 

(c) J r go(x, y)dfio(x) = for each y ET. 

(d) go(x, x) is a constant independent of x. 

Recall also that fx (with no subscript) is the admissible measure for the divisor D = K . By 
[21 Thm. 14.1] and the discussion in § |2.2[ we have the following further properties of /io and 
9o(x,y): 

go(x,y) = —-r(x,y) + C, for some constant C, (2.4) 



1 
1 



/i = - A, (r(x, y)) + 5 y , for any y 6 V, (2.5) 



fi = -(2/io + M. (2-6) 
2^ 



Now define a new function / : F — > R by 



1 f 1 

/(») = 2 y r ( x ^y) d ^(y) - ^g r ( x ' K )- 



We may rewrite / using (2.4) and (2.6) to get 

/(*) = ~y &o(x,y) (2 Mo (y) + 5 K {y)) +C- ^r{x,K) 

= -7r9o{x, K) + C- -7-r(x, K), by property (c). 
2# 4^ 

Taking the Laplacian of both sides with respect to x and applying property (b) and equa- 
tion (2.5) gives 

A,/ = (6 K - (2g - 2)/i ) - ]- (2(2g - 2)/i - 28 K ) = 0. 
2# 4g 

This implies that / is a constant, say / = A. Integrating / against \x gives 

A = ! f{x)dfi(x) = a (T) - — e (T) . 



2 The object /xq is called the canonical measure by Baker and Rumely p, §14]. 
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On the other hand, integrating against 5k yields 

(2g - 2) A = f(K) = le (F) - ±+(K,K). 
Solving these last two equations for a (T) completes the proof. □ 



Combining Proposition 2.3 with the definition of tp (r) immediately gives 



Corollary 2.4. For any pm-graph F of genus g, we have 

3. Reducing the Proof of Theorem 11.21 to Graph Theory 

In this section we give the argument that allows us to pass from a statement about algebraic 
curves to one about polarized metric graphs. While the argument is not difficult, it is new to 
the literature and relies on a recent theorem of Zhang. As such we give a detailed description. 
The reader may want to look at the beginning of §1.2| for references on semistable reduction 
and at [SI pp. 9-10] for further details on the theory of reduction graphs. 

3.1. Reduction Graphs. Suppose k is an algebraically closed field of characteristic zero 
and K is the function field of a smooth curve over k. Let Y be the proper smooth fc-curve 
with function field K. The points Y(k) are in bijective correspondence with places of K - 
i.e., discrete valuations of the function field K. For a point y G Y(k), let us write K y for the 
completion of K with respect to the corresponding discrete valuation. 

Let C be a proper smooth geometrically connected curve over K of genus g > 2. By the 
Semistable Reduction Theorem there exists a finite extension E/K y such that Ce = C x k E 
has semistable reduction. That is, if &e is the valuation ring of E, then there is a flat 
proper ^-scheme with generic fiber Ce and with semistable special fiber. Let us take 
/ : ^e — > Spec (Je to be the minimal regular (semistable) model of Ce- 

In what follows we identify y with the closed point of Spec <^£. Define a combinatorial 
graph G y with vertex set in bijection with the set of irreducible components of and 
with an edge between vertices p,p' for each node shared by the corresponding irreducible 
components Z, Z' . In particular, a vertex of G y admits a loop edge for each singularity of the 
corresponding irreducible component of f~ l (y). Endow G y with the structure of a weighted 
graph by defining £(e) = 1/[E : K y ] for each edge e. (As k is algebraically closed, [E : K y ] 
is the ramification index of the extension.) Define T y to be the metric graph induced by the 
weighted graph G y . Define q : T y — > Z as follows. Set q(p) = for all p £ V{G y ). For any 
vertex p of G y define q(p) to be the geometric genus of the irreducible component of f^iy) 
corresponding to the vertex p. Finally, set T y = (T y ,q). The key point of this construction 
is that, by virtue of the multiplicativity of ramification indices, the object T y is independent 
of the choice of extension E/K y . This follows for example from [51 Fact 2.5] J^] 

For example, if f~ l (y) is smooth, then Y y = {p} is a point, and q(p) = g, the genus of C. 
Conversely, if T y consists of a single point, then f~ l {y) is smooth. 



3 As the field of constants k is algebraically closed in our setting, the approach of Chinburg / Rumely in 
[5] can proceed without their notion of "well-adjusted" regular models. They require this extra technical 
tool in order to work with more general global fields. 
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Proposition 3.1. For each y e Y(k), the reduction graph T y = (T y ,q y ) is a pm-graph of 
genus g, where g is the genus of the curve C . 

Proof. To check that T y is a pm-graph, we only need to verify that its canonical divisor is 
effective. To avoide confusion, we will write D y for the canonical divisor of T y (the symbols 
K and K y being already in use). Recall that for p G T y the order of D y at p is given by 

ord p (A,) = 2q(p) + v(p) - 2, 

where v(p) is the valence of the point p. Let p±, . . . ,p n be the points of T y corresponding to 
irreducible components Z\, . . . , Z n of the fiber f~ 1 {y), numbered accordingly. It is clear from 
the definitions that ord p (D y ) = for any p ^ pj. The remaining cases can be calculated via 
intersection theory on the surface If ujf = uj^ e /^ e is the relative dualizing sheaf, then 
we claim 

ord Pi (A,) = Zi.Uf. (3.1) 

By the work of Arakelov, the sheaf Uf is numerically effective [81 Cor. 9.3.26]. So for any 
irreducible curve Zj C f^ 1 {y) we have Z^.Uf > 0, which shows the effectivity of D y . 

We now prove (3.1 ). Without loss of generality we may suppose i — 1, and we write p — p\ 
and Z — Z\. Also, let G y be the combinatorial graph dual to the fiber f^ 1 {y) as constructed 
above. We can compute the arithmetic genus of Z by the adjunction formula for regular 
fibered surfaces (8J Thm. 9.1.37]: 

*<z) = i + \ (z> + zm,) (3 2) 

=^ Z.oj f = 2p a (Z) - Z 2 - 2. 

The blowing-up formula relating arithmetic and geometric genus for a curve on a fibered 
surface behaves especially nicely in this case due to the presence of only nodal singularities. 
We have 

p a (Z) = p g {Z) + #{nodes of Z} 

— q(p) + #{loop edges of G y at p}. 

Define vi(p) (resp. v n {p)) to be the valence of G y at p contributed by loop edges (resp. 
non-loop edges). So v{p) = vi(p) + v n {p). Each loop contributes 2 to the valence at p, so 

Pa(Z) = q{p) + -Viip). (3.3) 

Next note that Z ~ — J2i>i %i since [f^ 1 (y)} = J2i ls a principle divisor. The intersec- 
tions of the Zi's with Z for % > 1 are in bijective correspondence with the non-loop edges of 
G y at p, so 

Z 2 = — Zi.Z = — #{non-loop edges at p} = —v n {p). (3.4) 

i>l 



Finally, we combine (3.2), (3.3), and (3.4) to find 

Zwj = 2 \q(p) + i«,(p)^ + u n (p) - 2 = 2g(p) + v(p) -2 = ord p (D y ), 
which is what we wanted. 
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To compute the genus of T y , we calculate the degree of its canonical divisor: 

2g (T y ) ~2=J2 ord p (D y ) = ^ord Pi (^) 

per y i 

i 

= deg{u CE ) =2g-2, 

where u>c E is the canonical sheaf on C^. Hence the genus of T y is g. □ 

3.2. The Work of Zhang. We continue to assume that C is a proper smooth connected 
curve of genus g > 2 over the function field K with field of constants k. Write C = C x K K . 
We begin with the following result of Zhang: 

Theorem 3.2 (Zhang, [T71 Thm. 5.6]). For any divisor D on C of degree 1, 

(uJ a ,uJ a ) ( T\ t ( K c 



a'(D) > v ~ a, ~ a ; + 1 - - )h [D - 



A(g-l) V 9 J V *9- 

where uj a is the admissible relative dualizing sheaf associated to C and Kc is a canonical 
divisor on C . 

We want to relate the admissible intersection number (u a , uj a ) to invariants of a global 
semistable model of C. Choose a finite extension K'/K such that Ck 1 = C K' has 
semistable reduction, and write d = [K f : K\. Let Y' be the proper smooth curve over k 
with function field K' . Define / : X' — > Y' to be the minimal regular (semistable) model of 
Ck 1 - Notice that we are working with a global semistable model of Ck> in contrast to the 
local model used in the previous section. Let (jJx'/y 1 be the relative dualizing sheaf. 

When / : X' — Y' is smooth and C is not a constant curve, it turns out that 



/ x ^X'/Y'-^X'/Y' . 12 
{UJ a , L0 a ) = > 



d d ' 

by a result of Parsin [12]. Then Zhang's result and positivity of the canonical height implies 

• t 'm^s i^a^a) . 3 

ml a (D) > — r > — r, 

DemvHc) 4(0-1) -d(g-iy 



which proves Theorem 1.2 when / is smooth. 

The rest of the discussion applies whether or not / is smooth, but it will be of greatest in- 
terest when / has some singular fibers. To relate the quantities a'(D) to our graph invariants 
in this case, we use the following recent result: 

Theorem 3.3 (Zhang, [181 Cor. 1.3.2 and §1-4]). With the notation above, 

y yeV(fe) 



Recall from £2.3 



that cp (T y ) = whenever T y consists of a single point. As C is smooth 
over K, there exists an affine subscheme U C Y and a proper smooth morphism X — ► U 
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with generic fiber C. Hence T y is a single point for all y G U(k), which means the above 
sum over points of Y is actually finite. 

Continuing with the notation above, let a : Y' — > Y be a morphism of proper smooth 
curves realizing the extension of function fields K' /K. For y' 6 F' (A;), we write r y / for the 
reduction graph of Ck> at the point y'. We need to relate if (IV J and <£> (r y ) whenever 
a(y') = y. In fact, we claim 



if (T y ,) = e y rcp (T y ) 



(3.5) 



where e y > is the ramification index of a at y'. The point is that T y > is defined relative to 
the base field K', so the lengths of its edges will differ from those of T y by exactly the 
ramification index. 

More precisely, let E = K' y , be the completion of K' with respect to the discrete valuation 
corresponding to y'. Then X' Xy< Specif is the minimal regular (semistable) model of 
Ce = Ck' x k'E because the fiber of X' over y' is unaffected by base change to the completed 
local ring Thus we can define F y i using this model. Notice that all of its segments will 
have length 1 because we did not need to make a finite extension of E = K' y , in order to 
obtain a semistable model over K'. On the other hand, K y C E since y' lies above y, and 
so we can define T y by passing to the extension E/K y . The degree [E : K y ] is equal to the 
ramification index e y i of a : Y' — > Y at y' . The underlying combinatorial graph structure of 
T y will be identical to that of T y > , but now a segm ent of T y will have length 1/[E : K y \ = l/e y . 

that if a pm-graph r' is obtained from T by 
= Xip (r). The claim (3.5) follows. 



4.6 



One can see from the proof of Proposition 
scaling all segments by the same quantity A, then \ T 

Next we need to know how the lengths £i (T y r) relate to the singular indices Si = 6i(X'/Y'). 
These notions are essentially dual to each other: a node of type % in the fiber f~ 1 (y') 
corresponds to a segment e of length 1 in T y >. Each point of the segment e (aside from the 
endpoints) is a point of T y > of type %. This correspondence allows one to verify that 

£i(?y>)=6 i (X'/Y'), 0<i<[g/2\. (3.6) 

y'eY'(k) 

We are now ready to reduce the main theorem to a statement about graph invariants. By 
positivity of the canonical height, we find 



inf a\D) > 

DeDiv^c) % 

> 



1) 



(Theorem 3.2 ) 



2(20 + 1 



2(20 + i; 



yGY(k) 

E \ E 

j/ey(fe) y'ey'W 

a(y')=y 



(Theorem 3.3) 



ej/ ^ (r v ) 



as 

a(y')=y 



(3.7) 



2d(2£+i; 



E 

?/'ey'( fc ) 



V ( r y') 



(by (3.5)). 



Thus we have reduced the problem about algebraic curves to a problem about pm-graphs. 
In the remainder of the paper, we will prove the following result: 
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Theorem 3.4. Let c(2) = i, c(3) 
0/ genus g = 2,3, or 4, we /iai>e 



81 ■ 



, and c(4) = ^ . T/ien for any polarized metric graph 



v (r) > c(g)£ (r) + £ 

ie(0,g/2] 



9 



U (r) . 



If we assume Theorem 3.4 for the moment, then the relations (3.7) and (|3.6|) show that 

2i(g - i) 



for a curve C/K of genus 2 < g < 4 

1 



inf a'(D) > - 
DGDiv^c) ~ 2d(2g + l 



1 



j/'ey'(fe) 



c(oKo(r y ,)+ £ 

ie (o, 9 /2] 



2d(2 5 + 1) 



ie(o, fl /2] 

Evidently this proves Theorem 1.2 Note that it gives a positive lower bound for a'(D) 
whenever there exists a singular fiber of / : X' — > Y 7 , and so the Geometric Bogomolov 
Conjecture is true for curves of genus at most 4 with bad reduction. 

4. Rapid Computation of ip (r) 

Let T = (r, q) be a pm-graph of genus g, and let us fix a model G. Enumerate the edges 
of G as ex, . . . , e m , and suppose these edges have lengths £1, . . . , £ m , respectively. By fixing 
the combinatorial type of G as well as the function q, we may view ip (r) as a function 
of £i,...,£ m . It is our goal now to further illuminate the nature of this function. We 
will use our new description to prove "continuity under edge contractions." This section is 
called rapid computation of ip (T) because we provide a reasonably efficient algorithm for its 



implementation in a computer algebra package. See Remarks 4.3 and 4.8 



Lemma 4.1. Let e be a segment of F with respect to the model G. Choose an orientation 
on G and use it to give an isometry y : [0,£(e)] e. For any vertex p of G, we have 

Y2 



r(p, y)dy 



6 V J 2 



r(p, e ) + r(p, e + ) 



Proof. By (2.5) and (2.6) in the proof of Proposition 2.3, we know that 

A y (r(p, y)) = 2gji + discrete masses. 



In particular, since we gave an explicit formula for /1 in {2.2 
- ^r(p,y(s)) = ) 



r(p,y(s)) 



F(e) 



-s 2 + As + B, 0<s< 



for some constants A and B independent of s. Substituting s = shows B = r(p,e ), and 
then substituting s = £(e) allows us to solve for A. Explicitly, we find 



r(p,y(s)) 



F(e) 



r(p, e + ) — r(p, e ) 



+ F{e) ) s + r(p, e~ 
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Now it is a simple matter of calculus to compute the integral in the lemma and arrive at the 
desired expression for it. □ 

Lemma 4.2. Let V = (T, q) be a pm-graph of genus g. Fix a model G of T with vertex set 
{pi, . . . ,p n } and edge set {ei, . . . ,e m }. Let K(pi) = v(pi) — 2 + 2q(pi) be the order of the 
canonical divisor at pi . Then we have 

* (r) = + ^ E F ^) 2 ^ 

9 id 69 k 



9 i k 

Proof. By definition, we have 

9 e P) = 9 J r{K, y)d^i{y) 

= 9^K(pi) I r(pi,y)dp:(y) 

i j \ j k k 

= Yq(Pj)K(p i )r(p ij p j ) + ^K(p i )^ ] -y^_ j r (pi,y)dy. 

i,j i k k Je h 

Inserting the formula from the previous lemma and simplifying completes the proof. □ 



Remark 4.3. Corollary 2.4 and Lemma 4.2 allow us to implement cp (T) in a computer 
algebra package. Indeed, all of the quantities involved are discrete in the sense that they 
depend only on the finite quantity of data contained in the weighted graph G and the 
function q. Computing the effective resistance is easy since it is essentially an inverse of the 
combinatorial Laplacian matrix. (Technically speaking, this is false since Q is singular.) We 
will push these ideas further momentarily in order to obtain more efficient algorithms for 
these computations. 

We continue with the notation from the beginning of this section. Define 

77(4,...,u = E (II'*)' ( 41 ) 

TCG \e k <£T J 

where the sum is over all spanning trees T of G. For example, r](l, . . . , 1) is the number of 
spanning trees of G. As each spanning tree is the complement of b± = dimnj i? 1 (r, M.) edges 
of G, we see that rj 6 ■ ■ ■ , £ m ] is homogeneous of degree b\. A useful alternative method 
for calculating rj is given by Kirchhoff's matrix-tree theorem. Let Q be the combinatorial 
Laplacian matrix for G. If the vertices of G are enumerated as p%, . . . ,p n , then Q is an n x n 
matrix whose entries are given by 

Se fe ={ Pl ,*}V lf*=j 
Qii= {-Y.e k = {Vl , Pj] Z~k l tii^j 



otherwise. 
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The first summation is over all non loop-edges e k containing the vertex pi and the second is 
over edges e k with vertices Pi and pj. 

The kernel of Q is generated by the vector [1,1,..., 1] T , and a simple consequence is that 
all of the first cofactors of Q are equal. (Use the easily verified fact Q ■ adj(Q) = 0.) Denote 
their common value by k*(G). Then the matrix-tree theorem says 

rn 

7 ] (£ l ,...,e m ) = K*(G)l[e k . (4.2) 

k=i 

Compare [4, §11.3, Thm. 12]. 

Lemma 4.4. Let T be a metric graph with first Betti number b\ = dimjR H 1 (r, M) ; and let 
G be a model of Y with vertex set {pi, . . . ,p n } and edge set {ex, . . . , e m }. For each pair of 
indices 1 < i,j < n, define 

Rij{£x, ■■ ■ ,£m) = r(pi,pj) rj(£x, . ■ .,£ m ). 
Then Rij is a homogeneous polynomial of degree b\ + 1 with integer coefficients. 

Proof. For i = j, we see immediately that = 0. By symmetry, it suffices to prove the 
result for fixed i and j with 1 < i < j < n. Set y = Pi and z = pj. By definition, 
r(y,z) = j z (y,y), where f(x) = j z (x,y) is the fundamental solution to the Laplace equation 
A x f = 5 y — 5 Z satisfying f(z) = 0. We may compute the value of j z (Pk,y) for any k using 
the discrete Laplace equation 

Qf = ei - e h fj = 0, (4.3) 

where is the A;th standard basis vector of 1R™ and f = J2 k fk e k- (Compare [U §5-6].) By 
definition, if f is the unique solution, then j z (Pk,y) = fk- 

The unique solution to the system of equations (4.3) gives rise to a solution of 

Q U) h = a, (4.4) 

where Q is the (n — 1) x (n — 1) matrix given by deleting the jth row and column from 
Q. The correspondence is given by 

j-l n 

f ^ h = f k u k + ^ fkU-k-i- 

k=l k=j+l 

Here is the kth standard basis vector of M n_1 . Since det(Q^) = k*(G) ^ by the 
matrix-tree theorem, the matrix Q^> is invertible and the solution to (4.4) is unique. 

Define to be the matrix with the ith column replaced by the zth standard basis 
vector Uj. The above argument and Cramer's rule shows that if h is the solution to (4.4), 

det(Q^) _ det(Qp } )nr=i4 _ det(Q?'>) UZi 4 



r(pi,Pj) = hi 



det(QO)) det(Q0'))nr=i4 v(h, 



R lj (£x,...,£ m ) = det(Q?)]l 



fc=i 



To complete the proof of the lemma, we must show that det(Q\ ■ ) YYk =1 4 is a homogeneous 
polynomial of degree bx + 1 in the lengths £x, ■ ■ ■ , £ m with integer coefficients. 
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To see that R^ is a polynomial, we assume that i = 1 and j = 2, perhaps after relabeling 
the vertices. Let G' be the graph given by fusing the vertices pi and p%. If P12 is the image 
of pi and p 2 in G', then we may write V{G') = {p±2,P3, • • • ,Pn} and E{G') = E{G). Let Q' 

be the Laplacian matrix of G' . As the edges adjacent to p 3 , . . . ,p n are unaffected by fusing 

(2) 

Pi and P2, the lower right (n — 2) x (n — 2) submatrixof Q agrees with that of Q\ . Hence 

(2) 

det(Q\ ) = k*(G ) in the notation preceding Lemma 



4.4 



It follows that 



%(4, . . . ,40 = det(QS 2) ) f]4 = ?7g'(4, • • • ,40: ( 4 - 5 ) 



fe=i 



where 77c is the polynomial associated to the graph G' as in (4.2). The definition (4.1) 
of t]g> shows that it is a polynomial in 4 •• • ,4i with integer coefficients. Moreover, the 
complement of any spanning tree in G' consists of bi + 1 edges, and so is homogeneous 
of degree b\ + 1. □ 

Lemma 4.5. Let V be a metric graph and G be a model of V with vertex set {pi, . . . ,p n } 
and edge set {e±, . . . , e m }. For an edge e k = {Pi,Pj} of G, set R k = Rij to be the polynomial 
defined in the previous lemma. Then £ k divides R k as polynomials in Z[4 . . . ,£ m ] and 

= r}(£i, .. .,£ m )\i k =o- 

Proof. The proof is similar to that of the previous lemma, except that we will contract the 
edge efc, rather than fusing its endpoints. After reordering the edges and vertices if necessary, 
we may assume that the edge of interest is e\ = {pi,P2}- Using the notation and strategy of 
the previous proof, we find 

Rliiu -- im) =det(Q?)f[£ k . 

1 k=2 

Define a new graph Gje\ given by contracting the edge e\. Let p\2 be the image of the 
vertices p\ and p 2 in G/e%. Then we may make the identifications E(G/e\) = E(G) \ {ei} 
and V(G/ei) = {pi2,P3, ■ ■ ■ ,p n }- Let Q{G/ei) be the Laplacian matrix of G/e\. It is not 
difficult to see that the lower right (n — 2) x (n — 2) sub matrix of Q(G /ei) agrees with that 



of Q[ , and so det(Qi ) = n*(G/ei) in the notation preceding Lemma 4.4 In particular, 

Rl(£l, ..-,40 /„ n V 

- B = VG/aih, •• • ,£m), 



where r]G/ ei is the polynomial associated to the graph Gje\ as in (4.2). 

To complete the proof, we must show that r] G / ei (£ 2 , . . . , £ m ) = r/(0, £2, ■ ■ ■ , 4)- But this 
follows from its definition and the bijective correspondence between spanning trees of Gje\ 
and spanning trees of G containing the edge ej. □ 

Proposition 4.6. Let T = (T, q) be a pm-graph of genus g > 1 and first Betti number 
b\. Fix a model G of V . With the notation above, there exists a homogeneous polynomial 
uj\ G Q[£i, ■ ■ ■ ,£ m ] °f degree 2b\ + 1, depending only on the combinatorial type of G and the 
function q, such that 

9-^ tt f^\ ^l(4, •••,4) 



V ( r ) = -1—^ V ) ~ 



6g ' 77(4, •••,4) ; 
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Proof. Write 77 = rj(£i, . . . ,£ m ) for simplicity. Let us begin by defining uj\ as the following 
function of £ x , . . . , £ m : 



1 ■ ■ ■ 1 ^mj 



65 



v 2 -v (r) v 2 - 



We must prove that uj\ is a homogeneous polynomial of degree 2b\ + 1 with rational coeffi- 
cients. 

As £{Y) = £1 + ■ ■ ■ + £ m , it is clear that £(T) rj 2 is homogeneous of the corr ect degree, so 
we are reduced to showing that <p (r) rj 2 has the same property. By Corollary 



2.4 



it suffices 



to verify that e (r) rj 2 and r(K,K) rj 2 are homogeneous degree-(2&! + 1) polynomials with 



rational coefficients. The latter is true by Lemma 4.4 

By Lemma 4.2 and the fact that F(ek) = 1 — r(ek)/£k, we have 



(r) rf^T, * (ft)* (Pi)*« + S^E^-?V^ 
9 i,j 69 k v tk ' 



Referring once more to Lemmas AA and 4J3, we see that this last expression is indeed a 
polynomial of the correct degree in the lengths £%, . . . , £ m with rational coefficients. □ 

Remark 4.7. Having computed a number of examples, it seems to be the case that the 
polynomial U\ defined in the previous proposition is divisible by the polynomial rj. While we 
cannot prove this in general, a strategy to prove it in the special case of hyperelliptic graphs 
is suggested by [TBI §3.3]. 

Remark 4.8. The results of this section and their proofs indicate the algorithms we have 
implemented in Mathematica. In all of these algorithms, the input is a weighted graph. If 
for example one wanted to calculate ip (T) for a metric graph, the best way to do it is to pick 
a model of T with very few vertices, and then apply these implementations to the model. 

1. The algorithm used to calculate effective resistance is given by the beginning of the 



4.6 



proof of Lemma |4.4 

To calculate e (r), we use the formula given by the proof of Proposition 
To calculate tp (r) , we use the proof of Proposition 4.6 to compute the polynomial uj\ 
and then apply the formula in the statement of the proposition. See the next remark 
for an explicit formula for uj\. 



Remark 4.9. For easy reference and later use, we now give an explicit formula for the poly- 
nomial associated to a metric graph V = (T, q) of genus g. This formula is easily derived 
from the arguments in this section. For simplicity, we use slightly different notation here; 
let R(x,y) = rj ■ r(x,y). Also, we write e = {e~,e + } for some choice of orientation on each 
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edge e. Then 



5^7-2 
+ 



V 



R(e 



k > e k 



V 



Sat a - 1) K (Pi) R (Pn Pi) {*9K{Pj) ~ 2(5(? - 2)q(Pj 



R(Pi, e k ; 



8^(5-1) 

By symmetry, this formula is independent of the choice of edge orientations. 

5. The Proof of Theorem 13.41 



Our goal for this section is to give the proof of Theorem |3.4[ which we restate for the 
reader's convenience. 

Theorem. Let c(2) = ^, c(3) = and c(4) = ^. Then for any polarized metric graph of 
genus g = 2,3, or 4, we have 

2i(g-i) 



v (r) > c(g)£ (r) + J2 

ie(o, s /2] 



9 



-ti (r) . 



Remark 5.1. The constant c(2) is sharp, but c(3) and c(4) are not optimal due to the use of 
a wasteful estimate in Proposition 5.16 below. For example, further computational evidence 
suggests that one could instead use c(3) = 17/288, but this has not been rigorously verified. 

The inequality in the theorem is a special case of the following more general conjecture: 

Conjecture 5.2 ([181 Conj. 4.1.1]). Let g > 2 be an integer. There exists a positive constant 
c(g) such that for any pm-graph T of genus g, 

2%{g 



<p (r) > c( 9 )e (r) + £ 

ie(o, s /2] 



9 



-ti (r) . 



Remark 5.3. Based on empirical evidence, we assert further that the conjecture should hold 
for polarized metric graphs of arbitrary genus g > 2 with 

9-1 



c(g) 



27g 



This is discussed at the end of Section 15.31 

Our plan is to show that, in order to prove the conjecture, it suffices to prove it for a 
finite simple class of graphs of genus g. (See Proposition 5.9 ) This will require a number 
of reduction steps. In the first part, we give a number of results on the continuity of certain 
metric graph invariants with respect to edge contraction. In the second part we show that, in 
order to give a lower bound for cp (T) for all pm-graphs of a given genus g, it suffices to give 
a lower bound when r = (r, q) is an irreducible cubic pm-graph of genus g with q identically 
zero. The basic idea is to reduce first to the case where V has no points of type i > by 
an additivity result on its components. Then we reduce to the case q = by adjoining 
circles along the support of q. Finally, we observe that an arbitrary metric graph with these 
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properties can be obtained from an irreducible cubic one by letting some of the edge lengths 
tend to zero. This last step uses the continuity results in the first part. A similar strategy 
was used in [HI H5J US] to reduce their respective questions to the case of cubic graphs. In 
the third part we give the strategy used to complete the proof of Theorem ^4 for graphs of 
a fixed genus. It requires a bit of difficult computation which we can accomplish when the 
genus is at most 4. In the fourth part we summarize these calculations, having performed 
them in Mathematica. The Mathematica code and the actual calculation notebooks close 
out the article. 



5.1. Continuity with Respect to Edge Contraction. In this section we restrict our 
attention to pm-graphs r = (T, q) with q = and genus g > 2. For simplicity, we will write 
T = (r, 0). The results hold for an arbitrary pm-graph with only slight modification, but we 
have no need for the general case. 

Let T 1 = (r 1; 0) be a pm-graph, and let G be a model of IY Fix an edge e\ of G that is 
not a loop edge, and let us identify it with the corresponding segment of F\. Define a family 
of pm-graphs Yt — (J^ti 0) as follows. For each t > 0, the combinatorial structure of is 
given by G, but the length of the segment e\ is t. For t — 0, let To be the metric graph given 
by contracting the segment e\ to a point; the weighted graph G/e\ serves as a model of r . 
Let ei = {^1,^2}, and let p±2 be the image of the segment e\ under the contraction Y\ — > T . 
As ei is not a loop edge, the family T t defines a homotopy from Ti to r . In particular, g(T t ) 
is constant for all t. 

Proposition 5.4. With notation as in the previous paragraph, 

V (T ) = Jim^ (T t ) . 

We will spend the rest of the section proving this result. Our strategy will be to use 
Proposition 4.6| and prove each of the quantities that appears there behaves well as t tends 



to zero. Each of these will be a separate lemma. 
Lemma 5.5. 

Kml(T t ) = £(r ). 

Proof. The definitions immediately give £(T t ) = t + £(Tq). □ 

Suppose the vertex and edge sets of G are given by {p 1 , . . . ,p n } and {ei, . . . , e m }, re- 
spectively. Let l\,ii,...,l m be the segment lengths of T t , where t\ = t. We also have 
E{G/e x ) = E(G) \ {ej and V(G/ ei ) = {p 12 ,p 3 , . . . ,p n }. Define 



EtcgIL^t^ if * > 0, 



where the summations are over spanning trees T. Note that ^(r t ) G Z[t,£ 2 , 
Lemma 5.6. 

jy( r *) = vFo) 
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Proof. Since the spanning trees of Gje\ are in bijection with the spanning trees of G con- 
taining the edge e 1; we see that as t — > 0, 

^ r *)= e n** +t £ 11** 

eiCTcG e fe (2T TcG e fe £T 

= e n^+°( i ) 

TCG/ei e^T 

= r ? (r ) + o(l). 

□ 

Let r t (x,y) denote the effective resistance between points x,y G IV Define 

R t (x,y) = Tj(T t ) ■ r t (x,y). 



We have already seen in Lemma |4T4 that R t (p,p') is a polynomial in the lengths t,£ 2 , ■ ■ ■ , 
for any vertices p, p' G G if t > (respectively p, p' G G/ei if t = 0). 



Lemma 5.7. 



lim RtiPiiPj] 



»/{i,j}c{l,2} 

Ro(pi2,Pj) ifi G {1,2} and j > 3 

Ro(Pi,Pn) if j G {1,2} audi > 3 

^RoiPiiPj) ifi,j>3. 



Proof. This result is intuitively obvious from the vantage point of circuit theory since the 
effective resistance behaves continuously under edge contractions. We will give a rigorous 



4.4 



(5.1) 



proof using the ideas from Lemma 4.4 

If pi = pj, then R t (pi,Pj) = for all t > 0. If i, j > 3, then R Q (pi,pj) = as well. Since 
we are in either the first or fourth cases of the lemma, we have proved the statement. 

For the remainder of the proof we assume pi ^ pj. By (4.5) in the proof of Lemma 
we find 

Rt(Pi,Pj) = VG'(tJ 2 , ■ ■ ■ Jm), 

where G' is the graph given by fusing the vertices Pi and pj. Thus 

limi?^,^) = 77 G /(0,£ 2 , ■■■Am) 

= e n 4 1 

TcG' e k <£T 

= e nt 

eiCTcG' e fe! z:T 

The first summation is over all spanning trees of G' while the second is only over the spanning 
trees containing e±. 

Let us now assume that {i,j} = {1,2}. Then the edge t\ becomes a loop edge in G', 
and so no spanning tree contains e\. The last expression in (5.1) must be zero, which is the 
desired result. 

Finally, we assume we are in one of the remaining three cases of the lemma. Then e\ is 
not a loop edge in G'. Let G'/ei be the graph derived from G' by collapsing the edge e\. 
It is naturally identified with the graph (G/ei)' given by first collapsing e\ and then fusing 
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the images of pi and pj. The spanning trees of G' containing e± are in bijection with the 
spanning trees of G'/ e i = (C/ e i)'- From ( |5.1[ ) we see that 



lim R t (p hPj ) = J2 H £ k 

erCTcG' e k <£T 

= e n«> 

TC(G/ei)' e^T 

= ^(G/ei)' (4, • • • , O (definition of r](G/ ei y) 
= Ro(pi,Pj), (by (|45[ again), 

where pi and are the images of Pi and pj, respectively, in the limit graph r . This completes 
the proof of the lemma. □ 



Let Ui(T t ) £ Q[t, £2, ■ ■ ■ 7 £ m ] be the polynomial defined in Proposition 4.6 with respect to 
the model G when t > (respectively Gje\ when t = 0). 



Lemma 5.8. In the setting above, 



limu;i(r t ) = wi(r ) 



Proof. Re mark 4.9 gives a formula for uji. Let K t be the canonical divisor for the graph Y t . 
As Lemma 4.5 and its proof imply r](T ) = Rt(e{ , e f)/t, we have 



wi(r t ) 



5g-2 
12g 



£(T t ) r)(T t f - t{rj(T t ) - n(r )) 2 - £ L(T t ) 



k>2 



R t(e k ,e+ 
4 



+ 



3r?(r-, 
807- 



1 -j 



5g-2 



(5.2) 



8g(g - 1) 

55-2 ri7/ ^ r /_,^ R t (e~,e k 



8g(g - 1) 



k>2 



4 



Rt(Pi,e k ) + Rtipuel 



We will treat each of these four expressions separately using the previous several lemmas. 
For the first expression of (5.2), we have 

hg-2 



I2g 



£(r t ) V (r t y-t (^-^(ro)) -E h( r 

fc>2 ^ 



5g-2 

I2g 



£(r ) v (r f-J2{v(r 



k>2 



+ o(l), 



as t — > 0. Thus we obtain the correct first expression for a;i(ro). Moving to the second 
expression, we note that K t {pi) = K (pi) for i > 3 and that 

K t ( Pl ) + K t (p 2 ) = v( Pl ) + v(p 2 ) - 4 = v(p 12 ) -2 = K (p 12 ). 
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Ignoring the leading factor on the second expression and using Lemma 5.7[ we see 
Y,K t { Pi )K t {p 3 )R t { Ph p 3 ) = 2K t (p 1 )K t (p 2 )R t (p 1 ,p 2 ) + £ KMKtipjRtfapj) 

+ 2 Kt(pi)K t (pi)Rt(pi,Pi) + 2 K t (p 2 )K t (pi)Rt(p2,Pi) 

i>3 i>3 

= K (pi)K (p j )R (p i ,p j )+2j2 [Ktfa) +K t (p 2 ^K (p i )R (p 12 ,p i ) +o(l) 

i,j>3 i>3 

= K {p i )K {p j )R {p l ,p j ) + 2Y, K o{Pi2)K { Pl )R {pi 2 ,Pi) + o{l). 

i,j>3 i>3 

Thus we find the second expression tends to the correct limit as t — > 0. 

The third expression in (5.2) vanishes as t tends to zero because of the factor T](T t ) — 
tj(Tq). The fourth expression tends to the correct limit using the same idea as in the second 
expression (i.e., group the terms corresponding to i = 1,2 separately). □ 



Proof of Proposition 5.4 By Proposition |4.6| and the previous lemmas, we have 



<? (r t ) = — nr.) - ^ = -^-/(r„) - ^ + o(i) 
= ¥>(r„)+ (i) 

as t — > 0. Here we have used the fact that r/(ro) ^ 0, which is evident by its definition and 
the fact that none of the other edge lengths is zero. □ 

5.2. Reduction to Irreducible Cubic PM-Graphs. In this section we show the following 
proposition: 

Proposition 5.9. Fix g > 2. Suppose there is a positive constant c(g) such that for all 
irreducible cubic pm-graphs Y = (T, 0) of genus g, 

y{T)>c{g)l{T). 
Then for an arbitrary pm-graph T, we have 

( p(T)> c (g)i (r)+ J2 2 ^ 1 ^(r)- 

*6 (0,9/2] 9 



In particular, it suffices to verify Conjecture 5J1 only for cubic irreducible pm-graphs T = 

(r,o). 

Our first step is to summarize some useful observations of Zhang [T8| §4.3]. Let T = 
(r, q) be an arbitrary pm-graph. Consider a finite collection of (closed, connected) metric 
subgraphs Tx,.. . ,T S such that T is the successive pointed sum of the IV For example, if 
T can be disconnected by removing a point p, then we may write r = Ti U T 2 for some 
subgraphs Tj such that I^ fl T 2 = {p}. In general, we may assign a pm-graph structure to 
each of the subgraphs Ti as follows. Let 7Tj : T — > Ti be the retraction map sending each 
point in T to the closest point in I\. Define a new function qi : Ti — > Z by 

q l {y) = dim R H 1 {n; 1 (y) J R) + £ q{ x ). 
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With these definitions, we have 



Lemma 5.10 ( |18[ Thm. 4.3.2]). Each pair Tj = (T^gj) is a pm-graph and g(T) = g(Ti). 
Moreover, 

^(r) = 5>( r *)- 

i 

Lemma 5.11 QlBJ Lem. 4.3.1, Prop. 4.4.1]). Let T = (T,q) be a pm-graph of genus g. Then 
T can be written as a successive pointed sum of finitely many subgraphs Ti and Ij such that 
each Ti is a maximal irreducible subgraph of T and each Ij is isometric to a closed interval 
for which q is identically zero on its interior. Moreover, 

*>(T) = I>(r.)+ £ ^-^.(r)- 

i ie(o,a/2] 



We extend our terminology slightly from {2.1 Let V = (T,q) be a pm-graph and p G T 



any point. Then p is of positive type if it is a point of type i for some i > 0. Equivalently, 
p is of positive type if p has valence 2, q(p) = 0, and T \ {p} is disconnected. Note that an 
irreducible metric graph has no points of positive type, but the converse is false. For example, 
the pointed sum of two circles can be disconnected by deleting their common point, but the 
type of this point is not defined. 

Lemma 5.12. Fix g > 2. Suppose there exists a positive constant c(g) such that for every 
pm-graph T = (T, q) of genus g with no points of positive type, 

<p(T)>c(g)e(r). 

Then for an arbitrary pm-graph of genus g, we have 

v(T)>c(g)l (T)+ 2 ^ 1 ^(T). 



ie(o, s /2] 



5.11 



we have 



Proof. Let T = (T, q) be an arbitrary pm-graph of genus g. By Lemma 

:i ie(o,ff/2] 

where Tj are the maximal irreducible subgraphs of T. By irreducibility, every point of Tj is of 
type or else the type is undefined. We may apply the hypothesis of the lemma to conclude 
that <f (Tj) > c(g)£(Tj). Evidently £ (T) = which completes the proof. □ 

Next we reduce to the case where q is identically zero. We first provide an explicit formula 
in the special case when T is a circle and q is supported at a single point. 

Lemma 5.13. Let C be a circle of length (-(C), let p G C be a point, and let g > 1 be an 
integer. Set C = (C, q) to be the pm-graph with q(p) = g — 1 and q(x) = for x ^ p. Then 

Proof. This computation is sufficiently simple to be done by hand and we leave it as an 
exercise for the reader. Alternatively, one could look at [HI Prop. 4.4.3], of which the 
present lemma is a special case. □ 
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Lemma 5.14. Fix g > 2. Suppose there exists a positive constant c(g) such that for every 
pm-graph F = (T, 0) of genus g with no points of positive type, 

cp(T)>c(g)£(T). (5.3) 

Then the above inequality holds for an arbitrary pm-graph of genus g with no points of positive 
type. 

Proof. We proceed by induction on the quantity ^2 xer q(x)- If XLer l( x ) = 0> then q = 
and (5.3) holds by hypothesis. Let us suppose that we have proved (5.3) holds for all pm- 
graphs T = (r, q) whenever Y has no points of positive type and XLer l( x ) — N- Choose 
r = (r, q) to be a pm-graph with N + 1 = ^2 xeT q(x) and no points of positive type. Then 
there is a point p G F with q{p) > 0. For each t > 0, define a new pm-graph T t = (F t , q t ) as 
follows: 

• The metric graph F t is obtained from V by gluing a circle of length t to the point p. 
Label the circle Ct- 

• Set 

(q(x) if x g r \ {p} 
q(p) - 1 if x = p. 
if x G C t n M 

Then F t is a pm-graph of genus g and XLer t l( x ) = Moreover, every point of C t \ {p} 
has type 0, and so T t has no points of positive type. By the induction hypothesis, the 
inequality (5.3) holds for Ft- 

Define C t = {C t , q') by setting q'(p) = g — 1 and q'(x) = otherwise. Then the induced pm- 
graph structures on the metric subgraphs F and C t of F t are given by F and C t , respectively. 
By the additivity described in Lemma 5.10[ we have 

<p (F t ) =<p(F)+<p (C t ) 

^ l p(F)= v (F t ) - <p (C t ) 

> c(g)£(Ft) — ^ — t by Lemma 5.13 

= c(g)£(F) + (c(g) - 9 -^j t. 

Let t — > to obtain the desired bound for cp (F) . This completes the induction step. □ 

Lemma 5.15. Fix g > 2. Suppose there exists a positive constant c(g) such that for every 
irreducible cubic pm-graph F = (F, 0) of genus g, 

¥>(T) >c(g)£(T). (5.4) 

Then the above inequality holds for an arbitrary pm-graph F = (F, 0) of genus g with no 
points of positive type. 

Proof. Let F = (F, 0) be a pm-graph of genus g with no points of positive type. It follows 
that if T \ {p} is disconnected for some p G F, then the valence of p is at least 3. If moreover 
T is cubic, then it must be irreducible. Indeed, if F \ {p} is disconnected, then p is a point 
of valence 3. There is a component of F \ {p} that contributes a single tangent direction 
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at p since three tangent directions must be shared among at least two components. Any 
point near p in this distinguished tangent direction will be a point of positive type, which is 
a contradiction. 

Now we proceed by reverse induction on the number of points of V of valence three, denoted 
N 3 (T). Observe that 

1 < Z>(p) - 2 ) 

per per 

u(p)=3 



N 3 (T) 



degK = 2g~2, 



and equality is achieved precisely when T is cubic. So if iV 3 (r) = 2g — 2, then (5.4) holds by 
the main hypothesis of the lemma. 

We now suppose that N 3 (T) < 2g — 2. There must exist a point p G T of valence v > 3. 
For each t > 0, we define a new graph I\ = (T t , 0) as follows. Delete the point po from T and 
compactify the resulting graph by adjoining new limit points X\, . . . , x v at the ends of each 
of the edges that previously met at p. If V \ {p } is disconnected, we further assume that x\ 
and X2 belong to distinct components of the compactification (after relabeling if necessary). 
The metric graph T t is obtained from this compactification by gluing x\ and X2 to one end 
of a closed segment of length t, and by gluing x 3 , . . . ,x u to the other end of the segment. 
(See Figure [2]) We call the segment eo = {pi,p 2 }, where pi is the image of x\ and 22, while 
P2 is the image of x 3 , . . . , x v . By construction Y t is connected and has no points of positive 
type for each t > 0. There is a canonical deformation retraction r t : T t — > T mapping the 
closed segment e to p - 

By construction the valence of p\ is 3. In particular, A^ 3 (Ff) > N 3 (T), so the induction 
hypothesis implies 

V {T t )>c(g)£(r t ) = c(g)£(T) + c(g)t. 
Therefore, to complete the proof it suffices to show that 

</?(r) =iim</?(r t ) . 



But this is precisely the content of Proposition 5.4| We have now completed the induction 
step. □ 



OQ 




CO. 



Figure 2. Here is an illustration of the construction of T t given in the proof 
of Lemma 5.15 in the case where the valence of p^ is v = 4. 



Proof of Proposition 5.S\ Apply the sequence of Lemmas 5.12 , 5.14 , and 5.15 



□ 



5.3. Strategy for the Proof of Theorem |3.4 

prove Theorem 



3.4 



We have seen in Proposition 5.9 that, to 
it suffices to prove it when T = (T, 0) is an irreducible cubic pm-graph 
of genus g > 2. Observe that such a graph admits a canonical model G° with precisely 
2g — 2 vertices. Indeed, we may take the set of valence 3 points as a vertex set, and the 
calculation in the proof of Lemma 5.15 shows this set has cardinality 2g — 2. This model 
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may have multiple edges, but it cannot have loop edges by irreducibility. Note further that 
g = #E{G°) - + 1, which shows 



#V(G°)=2g-2 #E(G°) = 3g-3. 

The next step is to list all possible combinatorial graphs that can arise as the canonical 
model G° as above, up to combinatorial isomorphism. One simple way to do it is to note 
that such a graph admits a maximal spanning tree consisting of 2g — 2 vertices and 2g — 3 
edges. It is easy to write down all of the trees of this type for small values of g. Given any 
such tree we build candidates for G° by adjoining g edges to it. Many of these will yield 
reducible graphs or graphs that fail to be 3-regular. Of the remaining candidates, many will 
be isomorphic. Rather than go through the details, we simply give the complete list of such 
graphs up to combinatorial isomorphism in Figure [3] below. 

Now let us suppose we have a combinatorial graph G° constructed as in the last paragraph. 



It is the canonical model of some metric graph T with edge lengths 4, • • 
3g — 3. By Proposition 4.6 , we know that 



where m 



9 



6g 



£(T) 



^i(4, 



e ) 2 



(5.5) 



Finally, we are reduced to the problem of bounding the quotient wi/rj 2 . The next propo- 
sition illustrates the strategy we will take for obtaining such a bound. 

Proposition 5.16. Let F = (T, 0) be an irreducible cubic pm-graph of genus g > 2. Let G° 
be its canonical model with edges ei, . . . , e m . Let a 2 o,nd a 3 be the second and third symmetric 
polynomials in the lengths 4, ■ • • , 4n- If there exists a positive constant A such that 



< A 



<?3(£l 



then 



<p(r)> 



9-1 



a 2 (ii, •••,40 
-A 39 ~t £(T). 



(5.6) 



6g ' 9(s - 1) 

Before beginning the proof, let us indicate how one might use the proposition to find the 



desired constant A. Clearing denominators in (5.6), we find the inequality 

A7 3 (4, • ■ ■ , 4nM4, • • • , tmf ~ ^l(4, • • • , 4nV 2 (4, • • • , 4n) > 0. 

The left side is a polynomial Yl a fa{A)£ a , where a is a multi-index and f a is a linear polyno- 
mial in A with rational coefficients. If there exists a value of A such that all of these linear 
polynomials are nonnegative, then of course f a (A)£ a > and the desired inequality 
holds. In order to find A, we solve the system of linear inequalities f a (A) > 0. As there can 
be thousands of such inequalities, we utilize Mathematica for this step. The computations 
are summarized in the next section, and the Mathematica notebooks are recreated afterward. 



Proof of Proposition 5.16 We begin by showing the inequality 



°A < m 

0~2 ~ 



3m 
27 



-£(T). 



(5.7) 



We use the following well-known generalization of the arithmetic-geometric mean inequality: 

l/m 





1 




1/2 




1/3 






> 




> 




> • • • > 


Cm 


Iff). 




L(s). 


n 

_ Vm/ _ 



where <7j is the ith symmetric polynomial in £1, 
find 

1/3 r 

G\ cr 2 



(T) (2) 



> 



c 3 

(3) 



, £ m . Using the first two inequalities, we 
1 2/3 



c 3 

(T) 



o- 3 



and (5.7) shows 



Rearranging and expanding out the binomial coefficients gives ( |5.7[ ), using <ji = ^(T). 
Now recall that the number of edges of G° is m = 3g — 3. Combining this fact with (5.6) 

• • • ,C) 



< A 



a 2 (£ 



4 < a^A^H- 

- 9(0-1) 1 J 



The proof is complete upon applying this last inequality to (5.5). 



□ 



We have now given the complete strategy for proving Theorem |3.4[ and all that remains 
is to provide the computations. To conclude this part, we give a general conjecture on the 
constant A in the proposition, and consequently a strong conjectural lower bound for (r) . 



Conjecture 5.17. With the notation of Proposition 5.16 , one may take 

7(9-1? 



A 



That is, Conjecture 5.2 holds with c(g) 



®g(3g - 5)' 
(g-l)/27g. 



This conjecture is born entirely from empirical data. It holds for all pm-graphs of genus 
g = 2,3,4, although it does not yield the smallest possible constant A in some cases. See 
Table ffl 

5.4 



of Thereom 3.4 



Computational Data. Here we provide all of the necessary data to complete the proof 
First, we enumerate the possible isomorphism classes of irreducible cubic 
combinatorial graphs G that can arise as the canonical model of a pm-graph T = (r, 0) of 
genus 2 < g < 4. As we have already indicated the strategy for finding all such graphs in 
the previous section, we simply list them in Figure [3j 

Table [T] summarizes the data obtained for the graphs in Figure [3j All of the computation 
was performed in Mathematica using the algorithms suggested in Section |4} The calculations 
are exact in the sense that no numerical methods were used. 

As there is only one graph of genus 2 in the table, we see immediately that Theorem |3.4 
holds with c(2) = 1/27. Both graphs of genus 3 yield the same lower bound for tp (r), so 
we may take c(3) = 2/81 in the theorem. For genus 4, two different lower bounds for ip (r) 
were obtained depending on whether the reductions in the previous section yield the graph 
r 4 2 . Of course we take the smaller of the two and conclude that c(4) = 1/36 will suffice in 
the theorem. 

Following Table [T] is the Mathematica code used for these computations. Some of the 
functions that appear there are useful for experimenting, but were not necessarily used to 
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produce the results in this article. Finally, we reproduce the Mathematica notebooks in 
which the computations were executed. 
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Genus 2 




Genus 3 




r 3 ,i 




Genus 4 





4,3 





Figure 3. Here we have a complete list of irreducible cubic pm-graphs with 
q = and genus at most 4, up to segment scaling and topological isomorphism. 
Equivalently, one can view this list as giving all possible combinatorial graphs 
(up to isomorphism) that can arise as the canonical model of such a pm-graph. 
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r 


A 


9-1 4 3g-5 
69 9(g-l) 


r 2 


5/12 


1/27 


r 3 i 


7/18 


2/81 


1 Q 9 


7/18 


2/81 


4,1 


3/8 


1 /36 Rj D D278 


j- 4,2 


7/26 


1 55 /28D8 D (1552 


r 4 , 3 


3/8 


1/36 


1^4,4 


3/8 


1/36 


^4,5 


3/8 


1/36 



Table 1 . The first column lists the names of the irreducible cubic pm-graphs 
given in Figure [3] The second column gives the smallest value of A that 
satisfies the inequality (5.6) of Proposition 5.16 while the third column gives 



the value of the constant provided by the conclusion of the proposition. 
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(* 
(* 
(* 



Metric Graph Calculation Tools 
Last update: 12.5.2008 
Written by X. W. C. Faber 



*) 



(* 



This collection of Mathematica functions is useful for computing many standard and new 
invariants of metric and pm-graphs. A brief description and usage of each function has been 
given below. Some of the detailed descriptions are left to the following two articles: 

[Zhang] S. Zhang. "Gross-Schoen Cycles and Dualising Sheaves." 
arXiv:0812.0371vl [math.NT], preprint, 2007. 

[Faber] X.W.C. Faber "The Geometric Bogomolov Conjecture for Small Genus Curves." 
arXiv:0803.0855v2 [math.NT], preprint, 2008. 

These functions were written for Mathematica 6.0.1.0. 
To load the functions, use the following commands: 

SetDirectory ["Current Local Directory"]; 
« "Combinatorica Mt 
« "GraphCalcTools - " 



ArithGenus::usage = " AritheGenus[g] returns the arithmetic genus of the graph g. It is the 
geometric genus plus the sum of all arithmetic weights on all of the vertices. This quantity 
is simply called the genus in [Faber] " 

AssignEdgeLengths::usage = " AssignLengths[g, L] returns a graph with the same 
underlying combinatorial structure as g, the same vertex weights, and with edge lengths given by 
the array L. This command makes it easier to define weighted graphs for computations using 
existing graphs " 

Assign Vertex Weights: :usage = " Assign Vertex Weightsfg, q] returns the graph with the same 
underlying combinatorial structure as g, the same edge weights, and with vertex weights given by 
the array q. This command makes it easier to define arithmetic graphs for computations using 
existing graphs " 

CanonicalDivisor::usage = " CanonicalDivisor[g] returns the canonical divisor of the graph g. It 
is given as a list {dl, dn} where di equals the valence of the ith vertex minus 2 plus twice the 
ith vertex weight. If g is not a polarized metric graph, then this returns the standard canonical 
divisor on a graph " 

EffRes::usage = " EffRes[g] gives the effective resistance matrix for the graph g. Its ij A th entry 
equals the effective resistance between vertices i and j of the weighted graph with n vertices" 

Epsilon::usage = " Epsilon[g] gives the epsilon constant defined in [Zhang] and [Faber] for the 
graph g. If g has unweighted edges it assigns weight 1 to them, and it assigns weight zero to any 
unweighted vertices. Note that this is the same as Moriwaki's or Yamaki's epsilon constant, 
although this is not obvious. See the references of [Faber] for further information on the work 
of Moriwaki and Yamaki" 

EqualEdgesLambda::usage = " EqualEdgesLambdafg] returns the Lambda constant for the graph 
g with equal edge lengths, total edge length 1, and no arithmetic weights. The constant Lambda 
was defined in [Zhang] " 



*) 
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EqualEdgesPhi::usage = " EqualEdgesfg] returns the Phi constant for the graph g with equal edge 
lengths, total edge length 1 , and no arithmetic weights " 

Eta::usage = " Eta[g] returns the Eta constant for the graph g. It is a homogeneous polynomial in 
the edge lengths of g of degree equal to the genus of g. See [Faber] for its definition " 

FosterCoeff::usage = " FosterCoeff[g] returns the m x 1 array of Foster coefficients for the graph 
g, where m is the number of edges. In the text of [Faber] the Foster Coefficient associated to 
edge e is denoted by F(e). It is given by the formula F(e) = 1 - r(e) / L(e). Here r(e) is the 
effective resistance between the endpoints of e and L(e) is the length of edge e " 

GeomGenus::usage = " GeomGenusfg] gives the geometric genus (first topological betti number, 
edges - vertices + 1, ... ) of the graph g " 

Lambda: :usage = " Lambda[g] gives the Lambda constant defined in [Zhang] for the graph g. If 
g has unweighted edges it assigns weight 1 to them, and it assigns weight zero to any unweighted 
vertices " 

Laplacian::usage = " Laplacian[g] returns the (vertex) Laplacian matrix of the weighted graph g " 

Omegal::usage = " Omegalfg] returns the polynomial omega_l defined in [Faber]. It is a 
homogeneous polynomial of degree 2g+l in the edge lengths of the graph " 

Omega2::usage = " Omega2[g] returns the polynomial omega_2, which is a homogeneous 
polynomial of degree 2g+l in the edge lengths of the graph, such that for a weighted graph g with 
arithmetic genus G, it satisfies the equality 

Lambdafg] = G / (8*G + 4) *TotalLength[g] + Omega2[g] / Eta[g] A 2. 
This quantity does not appear in either of the articles [Zhang] or [Faber] " 

Phi::usage = " Phi[g] gives the Phi constant defined in [Zhang] and [Faber] for the graph g. If g 
has unweighted edges it assigns weight 1 to them, and it assigns weight zero to any unweighted 
vertices " 

PMGraph::usage = " A graph g does not automatically have edge weights or vertex weights. 
PMGraph[g] represents the graph with the same edge and vertex set as g, but it assigns edge 
weight 1 to all edges that have no weights, and it assigns vertex weight zero to vertices that have 
no weight. PMGraph[g] = g if g already has weights assigned to all edges and vertices " 

Tau::usage = " Tau[g] gives the constant 'a' defined in [Faber] for the graph g. It is also the 
constant 'tau' defined in [Zhang]. For more details on how this constant relates to other graph 
invariants of Moriwaki and Yamaki, see the remarks in [Faber] following the definition of 'a'. If g 
has unweighted edges it assigns weight 1 to them, and it assigns weight zero to any unweighted 
vertices " 

TotalLength::usage = " TotalLength[g] returns the total length of the graph g. If g has unweighted 
edges, it assigns weight 1 to them. " 



*■) 



ArithGenus[g_Graph] := Module[ {Pg=PMGraph[g], m, n, q, i}, 

n = Dimensions[Pg[[2]]][[l]]; (* Number of Vertices *) 
m = Dimensions[Pg[[l]]][[l]]; (* Number of Edges *) 



33 



q = Sum[VertexWeight /. Pg[[2, i, 2]], {i,l,n}]; (* Total arithmetic weight *) 
Return[ m - n + 1 + q] 

] 



AssignEdgeLengths[g_Graph, L_] := Modulef { gNew=PMGraph[g], Len = L, m, i }, 
m = Dimensions[gNew[[l]]] [[1]]; (* Number of Edges *) 
Do[ 

gNew[[l,i, 2]] = (EdgeWeight -> 1/Len[[i]]), 
{i,l, m}]; 

ReturnfgNew] 

] 



Assign VertexWeights[g_Graph, q_] := Modulef { gNew=PMGraph[g], Weights=q, n, i }, 
n= Dimensions[gNew[[2]]][[l]]; (* Number of vertices *) 
Do[ 

gNew[[2,i, 2]] = ( VertexWeight -> Weights[[i]] ), 
{U,n}]; 

ReturnfgNew] 

] 



CanonicalDivisor[g_Graph] := Module[ {K, Pg=PMGraph[g], A,k, n, i, j }, 

n= Dimensions[Pg[[2]]][[l]]; (* Number of vertices *) 

K =Array[k, n ]; (* vector of multiplicities *) 

A = ToAdjacencyMatrix[Pg]; (* Adjacency Matrix *) 

Do[ 

K[[i]] = Sum[A[[i,j]],{j,l,n}] - 2 + 2* VertexWeight /. Pg[[2, i, 2]], 
{i,l,n}]; 

Return[K] 

] 



EffRes[g_Graph] := Module[{Pg = g, Q, ER, n, b, a, A, Al, i, j, k}, 

(* We need not call PMGraph because it is called inside Laplacian *) 
Q = Laplacian[Pg]; 

a = Det[Drop[Q, { 1 }, { 1 }]]; (* Weighted sum over spanning trees *) 

n= Dimensions[Pg[[2]]][[l]]; (* Number of vertices *) 

ER = Table[0, {i,n}, {j,n}]; (* Initialize Effective Resistance Matrix *) 

Do[ 

b=Table[0,{k,n-l}]; b[[i]]=l; 
A = Drop[Q, {j}, {j}]; (* Invertible Submatrix *) 
Al = Transpose[ReplacePart[Transpose[A], i->b]]; 
ER[[i,j]] = Simplify[ Det[Al] / a]; 
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ER[[j,i]] = ER[[i,j]];, 
{i, 1, n}, {j, i+1, n} ]; 

Return[ER] 

] 



Epsilon[g_Graph] := Module[{Pg=PMGraph[g], F, f, q, m, n, L, R, v, A, Eps, i, j, e}, 

n = Dimensions[Pg[[2]]][[l]]; (* Number of vertices *) 
m = Dimensions[Pg[[l]]][[l]]; (* Number of edges *) 

q = Table[VertexWeight /. Pg[[2, i, 2]] , {i, n}]; (* Array of arithmetic weights *) 

L = Table[ 1/EdgeWeight /. Pg[[l, i, 2]] , {i, m}]; (* Array of edge lengths *) 

R = EffRes[Pg]; (* Effective Resistance Matrix *) 

A = ToAdjacencyMatrix[Pg]; (* unweighted adjacency matrix *) 

v = Table[ Sum[A[[i,j]], {j,l,n}], {i,n}]; (* array of vertex valences *) 

(* Rather than call FosterCoeff, which requires another computation of the effective 
resistance matrix, we simply compute the Foster Coefficients here *) 
F = Array[f, m]; (* Foster coefficient array *) 
Do[ 

F[[i]] = 1 - R[[ Pg[[l,i,l,l]], Pg[[l,i,l,2]] ]]*EdgeWeight /. Pg[[l, i, 2]], 
{i,l,m}]; 

Eps = (l/ArithGenus[Pg]) * ( Sum[ q[[j]]*( v[[i]] - 2 + 2*q[[i]])*R[[ij]] , {i,l,n}, {j,l,n}] + 
Sum[ F[[e]]*( v[[i]] - 2 + 2*q[[i]])*( F[[e]]*L[[e]] / 6 + 
( R[[i, Pg[[l, e, 1, 2 ]] ]] + R[[i, Pg[[l, e, 1, 1 ]]]] ) 12 ),{i,l,n}, {e,l,m}] ); 

Return[Eps] 



EqualEdgesLambda[g_Graph] := Module[{Pg, m, L, i}, 

m = Dimensions[g[[l]]][[l]]; (* number of edges *) 
L = Table[l/m, {i,l,m}]; 
Pg = AssignEdgeLengths[g,L]; 
Return[Simplify[Lambda[Pg] ]] 



EqualEdgesPhi[g_Graph] := Module[{Pg, m, L, gNew=g, i}, 
m = Dimensions[gNew[[l]]][[l]]; (* number of edges *) 
L = Table[l/m, {i,l,m}]; 
Pg = AssignEdgeLengths[gNew,L]; 
Return[Simplify[Phi[Pg] ]] 



Eta[g_Graph] := Module[ {Pg = PMGraph[g], Q, ThisIsEta}, 

m = Dimensions[Pg[[l]]][[l]]; (* number of edges *) 
Q = Laplacian[Pg]; (* Laplacian *) 

L = Tablet 1/ EdgeWeight /.(Pg[[l, i, 2]] ), {i,l,m}]; (* Table of edge lengths *) 

P = Product[L[[e]], {e,l,m}]; (* Product of all of the edge lengths *) 

ThisIsEta = Expand[P*Det[Drop[Q, {!},{!}]]]; (* Weighted sum over spanning trees *) 
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Return[ThisIsEta] 

] 



FosterCoeff[g_Graph] := Module[ {Pg = PMGraphfg], m, R, F, f, i}, 

R = EffRes[Pg]; (* Effective Resistance Matrix *) 

m = Dimensions[Pg[[l]]][[l]]; (* Number of edges *) 

F = Array[f, m]; (* Foster coefficient array *) 

Do[ 

(* Foster coefficient of edge e is 1 - weight(e) x r(p,q), where p and q are the endpoints of e *) 
F[[i]] = 1 - R[[ Pg[[l,i,l,l]], Pg[[l,U,2]] ]]*EdgeWeight /. Pg[[l, i, 2]], 
{U,m}]; 

Return[F] 



GeomGenus[g_Graph] := Module[ {Pg=PMGraph[g], m, n}, 

n = Dimensions[g[[2]]][[l]]; (* Number of Vertices *) 
m = Dimensions[g[[l]]][[l]]; (* Number of Edges *) 
Return[ m - n + 1] 



Lambda[g_Graph] := Module[ {ThisIsLambda, Pg = PMGraph[g], G, Q, P, Eta, i, m}, 

G = ArithGenus[Pg]; 

Q = Laplacian[Pg]; (* Laplacian *) 

m = Dimensions[Pg[[l]]][[l]]; (* number of edges *) 

P = Product[l/(EdgeWeight /. Pg[[l,i,2]]), {i,l,m}]; (* Product of all of the edge lengths *) 
Eta = Expand[P*Det[Drop[Q, {1}, {1}]]]; (* Weighted sum over spanning trees *) 
ThisIsLambda = G / (8*G + 4) *TotalLength[Pg] + Omega2[Pg] / Eta A 2; 
Return[ThisIsLambda] 



Laplacian[g_Graph] := Module[{Pg = PMGraph[g], Diag, n, m, Q, i}, 

n= Dimensions[Pg[[2]]][[l]]; (* Number of vertices *) 

m = Dimensions[Pg[[l]]][[l]]; (* Number of Edges *) 

Q = Table[0, {i,n}, {j,n}]; (* Initialize Laplacian Matrix *) 

Do[ (* Cycle through the edges and update the matrix Q for each one *) 
If[ 

Pg[[l,i,l,l]] != Pg[[l,i,l,2]], (* Loop edges don't contribute *) 

Q[[Pg[[l,i,l,l]], Pg[[l,i,l,2]]]] = Q[[Pg[[Li,l,l]], Pg[[l,i,l,2]]]] - EdgeWeight/. Pg[[l,i,2]]; 
Q[[Pg[[l,i,l,2]], Pg[[l,i,l,l]]]] = Q[[Pg[[l,U,l]], Pg[[l,i,l,2]]]]; (* Symmetric Matrix *) 
Q[[Pg[[l,i,l,l]], Pg[[UU]]]] = Q[[Pg[[l,i,l,l]], Pg[[l,i,l,l]]]] + EdgeWeight /. Pg[[l,i,2]]; 
Q[[Pg[[l,i,l,2]], Pg[[l,i,l,2]]]] = Q[[Pg[[l,i,l,2]], Pg[[l,i,l,2]]]] + EdgeWeight /. Pg[[l,L2]]; 

], 

{U,m}]; 
Return[Q] 
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] 



Omega l[g_Graph] := Module[ {Pg = PMGraphfg], G, n, m, K, q, i, j, k, Q, P, Eta, L, ER, b, A, Al, e, 
ThisIsOmegal}, 

n = Dimensions[Pg[[2]]][[l]]; (* Number of vertices *) 
m = Dimensions[Pg[[l]]][[l]]; (* number of edges *) 
G = ArithGenus[Pg]; 

K = CanonicalDivisor[Pg]; (* Canonical Divisor *) 

q = Table[ VertexWeight /.(Pg[[2]][[i]][[2]] ), {i,l,n}]; (* Table of arithmetic weights *) 
Q = Laplacian[Pg]; (* Laplacian *) 

L = Tablet 1/ EdgeWeight /.(Pg[[l, i, 2]] ), {i,l,m}]; (* Table of edge lengths *) 

P = Product[L[[e]], {e,l,m}]; (* Product of all of the edge lengths *) 

Eta = Expand[P*Det[Drop[Q, {1}, {1}]]]; (* Weighted sum over spanning trees *) 

(* Now we compute the matrix ER: the effective resistance multiplied by Eta *) 
ER = Table[0, {i,n}, {j,n}]; (* Initialize Effective Resistance Matrix *) 

Do[ 

b=Table[0,{k,n-l}]; b[[i]]=l; 
A = Drop[Q, {j}, {j}]; (* Invertible Submatrix *) 
Al = Transpose[ReplacePart[Transpose[A], i -> b]]; 
ER[[ij]] = P*Det[Al]; 
ER[[j,i]] = ER[[ij]];, 
{i, 1, n}, {j, i+1, n} ]; 

ThisIsOmegal = (5*G - 2) / (12*G)*TotalLength[Pg]*Eta A 2 

+ Eta / (8*G*(G-1)) * Sum[ K[[i]]*ER[[i,j]]* (3*G*K[[j]] - 2*(5*G-2)*q[[j]]), {i,l,n}, {j,l,n}] 

- (5*G - 2) / (12*G)*Sum[ (Eta - ER[[Pg[[l, e, 1, 1 ]], Pg[[l, e, 1, 2 ]]]] / L[[e]] ) A 2*L[[e]], {e,l,: 

- (5*G - 2) / (8*G*(G-1)) * Sum[ K[[i]]*(Eta - ER[[Pg[[l, e, 1, 1 ]], Pg[[l, e, 1, 2 ]]]] / L[[e]] )* 

( ER[[i,Pg[[l, e, 1, 1 ]] ]] + ER[[i,Pg[[l, e, 1, 2 ]] ]] ), {i,l,n}, {e,l,m}]; 

Return[ThisIsOmegal ] 



Omega2[g_Graph] := Module[ {Pg = PMGraphfg], G, n, m, K, q, i, j, k, Q, P, Eta, L, ER, b, A, Al, e, 
ThisIsOmega2}, 

n = Dimensions[Pg[[2]]][[l]]; (* Number of vertices *) 
m = Dimensions[Pg[[l]]][[l]]; (* number of edges *) 
G = ArithGenus[Pg]; 

K = CanonicalDivisor[Pg]; (* Canonical Divisor *) 

q = Table[ VertexWeight /.(Pg[[2]][[i]][[2]] ), {i,l,n}]; (* Table of arithmetic weights *) 
Q = Laplacian[Pg]; (* Laplacian *) 

L = Tablet 1/ EdgeWeight /.(Pg[[l, i, 2]] ), {i,l,m}]; (* Table of edge lengths *) 

P = Product[L[[e]], {e,l,m}]; (* Product of all of the edge lengths *) 

Eta = Expand[P*Det[Drop[Q, {1}, {1}]]]; (* Weighted sum over spanning trees *) 

(* Now we compute the matrix ER: the effective resistance multiplied by Eta *) 
ER = Table[0, {i,n}, {j,n}]; (* Initialize Effective Resistance Matrix *) 

Do[ 

b=Table[0,{k,n-l}]; b[[i]]=l; 

A = Drop[Q, {j}, {j}]; (* Invertible Submatrix *) 
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Al = Transpose[ReplacePart[Transpose[A], i -> b]]; 
ER[[i,j]] = P*Det[Al]; 
ER[D,i]] = ER[[i,j]];, 
{i, 1, n}, {j, i+1, n} ]; 

Thisls0mega2 = (1-G) / (8*(2*G+l))*TotalLength[Pg]*Eta A 2 

+ Eta / (16*(2*G+1)) * Sum[ K[[i]]*ER[[i,j]]* (6*G*q[[j]] - K[[j]]), {i,l,n}, {j,l,n}] 
+ (G-l) / (8*(2*G+l))*Sum[ (Eta - ER[[Pg[[l, e, 1, 1 ]], Pg[[l, e, 1, 2 ]]]] / L[[e]] ) A 2*L[[e]], {e,l,m}] 
+ 3 / (16*(2*G+1)) * Sum[ K[[i]]*(Eta - ER[[Pg[[l, e, 1, 1 ]], Pg[[l, e, 1, 2 ]]]] / L[[e]] )* 
( ER[[i,Pg[[l, e, 1, 1 ]] ]] + ER[[i,Pg[[l, e, 1, 2 ]] ]] ), {i,l,n}, {e,l,m}]; 

Return[ThisIsOmega2] 



Phi[g_Graph] := Module[ {ThisIsPhi, Pg = PMGraph[g], G, Q, P, Eta, i, m}, 

G = ArithGenus[Pg]; 

Q = Laplacian[Pg]; (* Laplacian *) 

m = Dimensions[Pg[[l]]][[l]]; (* number of edges *) 

P = Product[l/(EdgeWeight /. Pg[[l,i,2]]), {i,l,m}]; (* Product of all of the edge lengths *) 
Eta = Expand[P*Det[Drop[Q, {1}, {1}]]]; (* Weighted sum over spanning trees *) 
ThisIsPhi = (G - 1) / (6*G) *TotalLength[Pg] - Omegal [Pg] / Eta A 2; 
Return[ThisIsPhi] 



PMGraph[g_Graph] := Module[{gNew=g, n, m, i}, 

n=Dimensions[gNew[[2]]] [[1]]; (* number of vertices *) 
m=Dimensions[gNew[[l]]] [[1]]; (* number of edges *) 

(* This loop adjusts the edge weights *) 
Do[ 

If[ 

Dimensions [gNew[ [ 1 ,i] ] , 1 ] [ [ 1 ] ]== 1 , 
gNew[[l,i]]= {gNew[[l,i,l]], Edge Weight ->1}; 

], 

{i,l, m}]; 

(* This loop adjusts the vertex weights *) 
Do[ 

Ifl 

Dimensions[gNew[[2,i]], 1][[1]]=1, 
gNew[[2,i]] = {gNew[[2,i,l]], VertexWeight ->0} 

], 

{i,l,n}]; 
ReturnfgNew] 



Tau[g_Graph] := Module[{Pg=PMGraph[g], F, q, m, n, L, R, v, A, T, i, j, e, f, k}, 

n = Dimensions[Pg[[2]]][[l]]; (* Number of vertices *) 
m = Dimensions[Pg[[ !]]][[!]]; (* Number of edges *) 
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q = Table[VertexWeight /. Pg[[2, i, 2]] , {i, n}]; (* Array of arithmetic weights *) 

L = Table[ 1/EdgeWeight /. Pg[[l, i, 2]] , {i, m}]; (* Array of edge lengths *) 

R = EffRes[Pg]; (* Effective Resistance Matrix *) 

A = ToAdjacencyMatrix[Pg]; (* unweighted adjacency matrix *) 

v = Table[ Sum[A[[i,j]], {j,l,n}], {i,n}]; (* array of vertex valences *) 

(* Rather than call FosterCoeff, which requires another computation of the effective 
resistance matrix, we simply compute the Foster Coefficients here *) 

F = Arrayfk, m]; (* Foster coefficient array *) 
Do[ 

F[[i]] = 1 - R[[ Pg[[l,i,l,l]], Pg[[l,i,l,2]] ]]*EdgeWeight /. Pg[[l, i, 2]], 
{U,m}]; 

T = (1 / ( 2*ArithGenus[Pg] A 2 ) )* ( Sum[ q[[i]]*q[[j]]*R[[i,j]], {U,n}, {j,l,n} ] + 

2* Sum[ q[[i]]* ( (l/6)*F[[e]] A 2 * L[[e]] + (l/2)*F[[e]]* 

( R[[i, Pg[[l, e, 1, 2 ]] ]] + R[[i, Pg[[l, e, 1, 1 ]]]] ) ), {e,l,m}, {i,l,n} ] + 

2*Sum[ F[[e]]*F[[f]]* ( (l/6)*(F[[e]]*L[[e]] +F[[f]]*L[[f]]) + 

(l/4)*( R[[Pg[[l, e, 1,2]], Pg[[l, f, 1, 2 ]] ]] + R[[Pg[[l, e, 1,1]], Pg[[l, f, 1, 1]] ]] 
+ R[[Pg[[l, e, 1, 2 ]] , Pg[[l, f, 1, 1 ]] ]] + R[[Pg[[l, e, 1,1]], Pg[[l, f, 1, 2 ]] ]] ) ), 
{e , 1, m}, {f, e+1, m}] + Sum[ L[[e]] / 6 * F[[e]] A 2 *(2 - F[[e]] ) ,{e,l,m} ] ); 

Return[T] 



TotalLength[g_Graph] := Module[{Pg=PMGraph[g], L, e, Total, m, i}, 

m = Dimensions[Pg[[l]]][[l]]; (* Number of edges *) 

L = Table[ 1/EdgeWeight /■ Pg[[l, i. 2]] , {i, m}]; (* Array of edge lengths *) 

Total = Sum[ L[[e]],{e,l,m}]; 

Return[Total] 



Print["Loaded ArithGenus, AssignEdgeLengths, AssignVertex Weights, CanonicalDivisor, EffRes, Epsilon, 
EqualEdgesLambda, EqualEdgesPhi, Eta, FosterCoeff, GeomGenus, Lambda, Laplacian, Omega 1, Omega2, Phi, 
PMGraph, Tau, TotalLength"] 

(* End File *) 
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Mathematics Notebook for the Graph r 2 



This notebook collects all of the calculations relevant to computing some constant A satisfying the hypothesis of 
Proposition 5.16 for the graph r=r 2 . 

In[2]:= SetDirectory [Current] ; 
<< "Combinatorica^ " ; 
<< "GraphCalcTools~ " ; 

Loaded ArithGenus, AssignEdgeLengths, AssignVertexWeights, CanonicalDivisor, 
EffRes, Epsilon, EqualEdgesLambda, EqualEdgesPhi, Eta, FosterCoeff, 
GeomGenus, Lambda, Laplacian, Omegal, 0mega2, Phi, PMGraph, Tau, TotalLength 

First, define the array of edge lengths and the combinatorial structure of the graph in question, and then define the 
relevant functions of the edge lengths. 

In[5]:= L = {a, b, c} ; 

g=Graph[{{{l, 2}}, {{1,2}}, {{1,2}}}, {{{-1,0}}, {{1,0}}}]; 

r = AssignEdgeLengths [g, L] ; 

ShowGraph[r] 



Out[8]= 



ln[9]:= r) = Expand [Eta [T] ] ; 

till = Expand [Omegal [r] ] ; 

o2 = SymmetricPolynomial [2 , L] ; 

a3 = SymmetricPolynomial [3 , L] ; 

The next step uses the strategy following the statement of Proposition 5.16 to create a list of linear inequalities in the 
variable A. The object PI is an array of the coefficients of the polynomial A*?7 A 2*cr3 - wl*<x2, viewed as a function of 
the edge lengths. Each coefficient is a linear polynomial in A. We verify that it is a 4x4x4 array. 

In[l3]:= PI = Coef f icientList [Expand [A * r) * 2 * a3 - a>l * a2] , L] ; 
Dimensions [PI] 

Out[14]= {4, 4, 4} 

The Do - loop creates a new list Coefflneq that consists of linear inequalities built from the nonzero entries of PI. The 
command "Reduce" solves this list of linear inequalities. 

In[15]:= Coef f Ineq = {}; 
Do[ 

Test = PI [ [a, b, c] ] ; 

If [Test + | | D[Test, A] + 0, Coef f Ineq = Append [Coef f Ineq, Test > 0] ] ; 
, {a, 4}, {b, 4}, {c, 4}] ; 
Reduce [Coefflneq] 

5 

Out[17]= A > 

12 
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Mathematics Notebook for the Graph r 3/1 

This notebook collects all of the calculations relevant to computing some constant A satisfying the hypothesis of 
Proposition 5.16 for the graph r=T3 i. 

First, define the array of edge lengths and the combinatorial structure of the graph in question, and then define the 
relevant functions of the edge lengths. 

In[18]:= L = Array [f , {6}]; 

r = AssignEdgeLengths [ 

Graph[{{{l, 2}}, {{1,3}}, {{1,4}}, {{2,3}}, {{2,4}}, {{3,4}}}, 
{{{0,0}}, {{0,1}}, {{-Sqrt[3] /2, -1/2}}, {{Sqrt[3]/2, -1/2}}}], L] ; 
ShowGraph [ 
r] 



Out[20]= 




ln[21]:= r) = Expand [Eta [T] ] ; 

col = Expand [ Omega 1 [r] ] ; 

o2 = SymmetricPolynomial [2 , L] ; 

a3 = SymmetricPolynomial [3 , L] ; 

The next step uses the strategy following the statement of Proposition 5.16 to create a list of linear inequalities in the 
variable A. The object PI is an array of the coefficients of the polynomial A*r] A 2*cr3 - oj1*<x2, viewed as a function of 
the edge lengths. Each coefficient is a linear polynomial in A. We verify that it is 6-dimensional array with 4 entries in 
each dimension. 

In[25]:= PI = Coef f icientList [Expand [A * r) " 2 * a3 - ul*a2], L] ; 
Dimensions [PI] 

Out[26]= {4, 4, 4, 4, 4, 4} 

The Do - loop creates a new list Coefflneq that consists of linear inequalities built from the nonzero entries of PI. The 
command "Reduce" solves this list of linear inequalities. 

In[27]:= Coef f Ineq = {}; 
Do[ 

Test = PI [ [a, b, c, d, e, f ] ] ; 

If [Test t | | D[Test, A] # 0, Coef f Ineq = Append [Coef f Ineq, Test t. 0] ] ; 
, {a, 4}, {b, 4}, {c, 4}, {d, 4}, {e, 4}, {f, 4}]; 
Reduce [Coefflneq] 

7 

Out[29]= A > 

18 
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Mathematics Notebook for the Graph T 3f2 

This notebook collects all of the calculations relevant to computing some constant A satisfying the hypothesis of 
Proposition 5.16 for the graph T=r^,2- 

First, define the array of edge lengths and the combinatorial structure of the graph in question, and then define the 
relevant functions of the edge lengths. 

In[30]:= L =Array[Y, {6}]; 

r = AssignEdgeLengths [ 

Graph[{{{l, 2}}, {{1,2}}, {{2,3}}, {{1,4}}, {{3,4}}, {{3,4}}}, 
{{{1,1/3}}, {{0,1/3}}, {{0,0}}, {{1,0}}}], L]; 
ShowGraph [ 
r] 



Out[32]= 




ln[33]:= rj = Expand [Eta [r] ] ; 

0)1 = Expand [ Omega 1 [r] ] ; 

o2 = SymmetricPolynomial [2 , L] ; 

a3 = SymmetricPolynomial [3 , L] ; 

The next step uses the strategy following the statement of Proposition 5.16 to create a list of linear inequalities in the 
variable A. The object PI is an array of the coefficients of the polynomial A*r] A 2*cr3 - oj1*<x2, viewed as a function of 
the edge lengths. Each coefficient is a linear polynomial in A. We verify that it is 6-dimensional array with 4 entries in 
each dimension. 

\n[37]:= PI = Coef f icientList [Expand [A * r) " 2 * a3 - ul*a2], L] ; 
Dimensions [PI] 

Out[38]= {4, 4, 4, 4, 4, 4} 

The Do - loop creates a new list Coefflneq that consists of linear inequalities built from the nonzero entries of PI. The 
command "Reduce" solves this list of linear inequalities. 

In[39]:= Coef f Ineq = { } ; 
Do[ 

Test = PI [ [a, b, c, d, e, f ] ] ; 

If [Test f | | D[Test, A] f 0, Coef f Ineq = Append [Coef f Ineq, Test > 0] ] ; 
, {a, 4}, {b, 4}, {c, 4}, {d, 4}, {e, 4}, {f, 4}]; 
Reduce [Coefflneq] 

7 

Out[41]= A > 

18 
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Mathematics Notebook for the Graph r 4/1 



This notebook collects all of the calculations relevant to computing some constant A satisfying the hypothesis of 
Proposition 5.16 for the graph T=T^i. 

First, define the array of edge lengths and the combinatorial structure of the graph in question, and then define the 
relevant functions of the edge lengths. 

In[42]:= L =Array[Y, {9}]; 

r = AssignEdgeLengths [ 

Graph[{{{l, 2}}, {{3,4}}, {{5, 6}}, {{2,4}}, {{4,6}}, {{2, 6}}, {{1,3}}, 
{{3,5}}, {{1,5}}}, {{{0,1}}, {{0,2}}, {{-Sqrt[3] /2, -1/2}}, 
{{-Sqrt[3], -1}}, {{Sqrt[3]/2, -1/2}}, {{Sqrt[3], -1}}}], L] ; 
ShowGraph [ 
r] 



Out[44]= 




ln[45]:= 1 = Expand [Eta [r] ] ; 

ul = Expand [ Omega 1 [r] ] ; 

a2 = SymmetricPolynomial [2 , L] ; 

a3 = SymmetricPolynomial [3, L] ; 

The next step uses the strategy following the statement of Proposition 5.16 to create a list of linear inequalities in the 
variable A. The object PI is an array of the coefficients of the polynomial A*r] A 2*cr3 - wl*cr2, viewed as a function of 
the edge lengths. Each coefficient is a linear polynomial in A. We verify that it is a 9-dimensional array with 4 entries in 
each dimension. 

In[49]:= PI = Coef f icientList [Expand [A * r) * 2 * a3 - ul*a2], L] ; 
Dimensions [PI] 

Out[50]= {4, 4, 4, 4, 4, 4, 4, 4, 4} 

The Do - loop creates a new list Coefflneq that consists of linear inequalities built from the nonzero entries of PI. The 
command "Reduce" solves this list of linear inequalities. 

In[51]:= Coef f Ineq = { } ; 
Do[ 

Test = PI [ [a, b, c, d, e, f, g, h, i] ] ; 

If [Test f | | D[Test, A] f 0, Coef f Ineq = Append [Coef f Ineq, Test > 0] ] ; 
, {a, 4}, {b, 4}, {c, 4}, {d, 4}, {e, 4}, {f, 4}, {g, 4}, {h, 4}, {i, 4}]; 
Reduce [Coefflneq] 

3 

Out[53]= A > - 
8 
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Mathematics Notebook for the Graph r 4 , 2 



This notebook collects all of the calculations relevant to computing some constant A satisfying the hypothesis of 
Proposition 5.16 for the graph T=T^2- 

First, define the array of edge lengths and the combinatorial structure of the graph in question, and then define the 
relevant functions of the edge lengths. 

In[54]:= L =Array[Y, {9}]; 

r = AssignEdgeLengths [ 

Graph[{{{l, 2}}, {{1,2}}, {{1, 6}}, {{2,3}}, {{3,4}}, {{3, 4}}, {{4,5}}, 
{{5,6}}, {{5,6}}}, {{{1,0}}, {{1/2, Sqrt[3] /2}}, {{-1/2, Sqrt[3]/2}}, 
{{-1,0}}, {{-1/2, -Sqrt[3] /2}}, { { 1 / 2 , -Sqrt [3] / 2} } } ] , L] ; 
ShowGraph [ 
r] 



Out[56]= 




ln[57]:= r) = Expand [Eta [T] ] ; 

oil = Expand [ Omega 1 [r] ] ; 

a2 = SymmetricPolynomial [2 , L] ; 

a3 = SymmetricPolynomial [3 , L] ; 

The next step uses the strategy following the statement of Proposition 5.16 to create a list of linear inequalities in the 
variable A. The object PI is an array of the coefficients of the polynomial A*?7 A 2*cr3 - wl*cr2, viewed as a function of 
the edge lengths. Each coefficient is a linear polynomial in A. We verify that it is a 9-dimensional array with 4 entries in 
each dimension. 

In[61]:= PI = Coef f icientList [Expand [A * r) * 2 * a3 - ail *a2], L] ; 
Dimensions [PI] 

Out[62]= {4, 4, 4, 4, 4, 4, 4, 4, 4} 

The Do - loop creates a new list Coefflneq that consists of linear inequalities built from the nonzero entries of PI. The 
command "Reduce" solves this list of linear inequalities. 

In[63]:= Coefflneq = { } ; 
Do[ 

Test = PI [ [a, b, c, d, e, f, g, h, i] ] ; 

If [Test f | | D[Test, A] f 0, Coef f Ineq = Append [Coef flneq. Test > 0] ] ; 
, {a, 4}, {b, 4}, {c, 4}, {d, 4}, {e, 4}, {f, 4}, {g, 4}, {h, 4}, {i, 4}]; 
Reduce [Coefflneq] 

7 

Out[65]= A > 
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Mathematics Notebook for the Graph r 4 , 3 

This notebook collects all of the calculations relevant to computing some constant A satisfying the hypothesis of 
Proposition 5.16 for the graph I^r^. 

First, define the array of edge lengths and the combinatorial structure of the graph in question, and then define the 
relevant functions of the edge lengths. 

In[66]:= L =Array[Y, {9}]; 

r = AssignEdgeLengths [ 

Graph[{{{l, 2}}, {{1,6}}, {{1, 6}}, 

{{2,3}}, {{2,5}}, {{3, 4}}, {{3, 4}}, {{4, 5}}, {{5,6}}}, 
{{{1,1}}, {{0,1}}, {{-1, 1}}, {{-1,0}}, {{0,0}}, {{1,0}}}], L]; 
ShowGraph [ 
r] 



Out[68]= 




ln[69]:= rj = Expand [Eta [r] ] ; 

0)1 = Expand [ Omega 1 [r] ] ; 

o2 = SymmetricPolynomial [2 , L] ; 

a3 = SymmetricPolynomial [3 , L] ; 

The next step uses the strategy following the statement of Proposition 5.16 to create a list of linear inequalities in the 
variable A. The object PI is an array of the coefficients of the polynomial A*77 A 2*cr3 - u>l*cr2, viewed as a function of 
the edge lengths. Each coefficient is a linear polynomial in A. We verify that it is a 9-dimensional array with 4 entries in 
each dimension. 

In[73]:= PI = Coef f icientList [Expand [A * r) " 2 * a3 - ul*a2], L] ; 
Dimensions [PI] 

Out[74]= {4, 4, 4, 4, 4, 4, 4, 4, 4} 



The Do - loop creates a new list Coefflneq that consists of linear inequalities built from the nonzero entries of PI. The 
command "Reduce" solves this list of linear inequalities. 

In[75]:= Coefflneq = { } ; 
Do[ 

Test = PI [ [a, b, c, d, e, f, g, h, i] ] ; 

If [Test f | | D[Test, A] f 0, Coef f Ineq = Append [Coef flneq. Test > 0] ] ; 
, {a, 4}, {b, 4}, {c, 4}, {d, 4}, {e, 4}, {f, 4}, {g, 4}, {h, 4}, {i, 4}]; 
Reduce [Coefflneq] 

3 

Out[77]= A > - 
8 
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Mathematics Notebook for the Graph r 4 , 4 

This notebook collects all of the calculations relevant to computing some constant A satisfying the hypothesis of 
Proposition 5.16 for the graph I^r^. 

First, define the array of edge lengths and the combinatorial structure of the graph in question, and then define the 
relevant functions of the edge lengths. 

In[78]:= L =Array[Y, {9}]; 

r = AssignEdgeLengths [ 

Graph[{{{l, 2}}, {{1,4}}, {{1, 5}}, 

{{2,3}}, {{2,6}}, {{3, 4}}, {{3, 5}}, {{4, 6}}, {{5,6}}}, 
{ {{-1, 1}}, {{0, 1}}, {{1, 1}}, {{2, 1}}, {{0, 0}}, {{1, 0}} }], L] ; 
ShowGraph [ 
r] 



Out[80]= 




For whatever reason, the command ShowGraph doesn 1 1 display the edge between the upper left and upper right vertices. 

In[81]:= 1 = Expand [Eta [T] ] ; 

0)1 = Expand [ Omega 1 [r] ] ; 

o2 = SymmetricPolynomial [2 , L] ; 

a3 = SymmetricPolynomial [3, L] ; 

The next step uses the strategy following the statement of Proposition 5.16 to create a list of linear inequalities in the 
variable A. The object PI is an array of the coefficients of the polynomial A*?7 A 2*cr3 - wl*cr2, viewed as a function of 
the edge lengths. Each coefficient is a linear polynomial in A. We verify that it is a 9-dimensional array with 4 entries in 
each dimension. 

In[85]:= PI = Coef f icientList [Expand [A * T) " 2 * a3 - a>l*a2], L] ; 
Dimensions [PI] 

Out[86]= {4, 4, 4, 4, 4, 4, 4, 4, 4} 

The Do - loop creates a new list Coefflneq that consists of linear inequalities built from the nonzero entries of PI. The 
command "Reduce" solves this list of linear inequalities. 

In[87]:= Coefflneq = { } ; 
Do[ 

Test = PI [ [a, b, c, d, e, f, g, h, i] ] ; 

If [Test t | | D[Test, A] # 0, Coef f Ineq = Append [Coef flneq. Test > 0] ] ; 
, {a, 4}, {b, 4}, {c, 4}, {d, 4}, {e, 4}, {f, 4}, {g, 4}, {h, 4}, {i, 4}]; 
Reduce [Coefflneq] 

3 

Out[89]= A > — 
8 
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Mathematics Notebook for the Graph r 4 , 5 



This notebook collects all of the calculations relevant to computing some constant A satisfying the hypothesis of 
Proposition 5.16 for the graph r=T4 5. 

First, define the array of edge lengths and the combinatorial structure of the graph in question, and then define the 
relevant functions of the edge lengths. 

In[90]:= L =Array[Y, {9}]; 

r = AssignEdgeLengths [ 

Graph[{{{l, 2}}, {{1,4}}, {{1, 5}}, 

{{2,3}}, {{2,4}}, {{3, 4}}, {{3, 6}}, {{5, 6}}, {{5,6}}}, 
{{{1,1}}, {{1/2,0}}, {{1,-1}}, {{3/2,0}}, {{-1,1}}, {{-1,-1}}}], L]; 
ShowGraph [ 
r] 

r 



Out[92]= 



I 

ln[93]:= r) = Expand [Eta [r] ] ; 

oil = Expand [ Omega 1 [r] ] ; 

o2 = SymmetricPolynomial [2 , L] ; 

a3 = SymmetricPolynomial [3 , L] ; 

The next step uses the strategy following the statement of Proposition 5.16 to create a list of linear inequalities in the 
variable A. The object PI is an array of the coefficients of the polynomial A*j] a 2*(t3> - wl*<x2, viewed as a function of 
the edge lengths. Each coefficient is a linear polynomial in A. We verify that it is a 9-dimensional array with 4 entries in 
each dimension. 

In[97]:= PI = Coef f icientList [Expand [A * r) * 2 * a3 - 0)1 *a2], L] ; 
Dimensions [PI] 

Out[98]= {4, 4, 4, 4, 4, 4, 4, 4, 4} 

The Do - loop creates a new list Coefflneq that consists of linear inequalities built from the nonzero entries of PI. The 
command "Reduce" solves this list of linear inequalities. 

In[99]:= Coef f Ineq = { } ; 
Do[ 

Test = PI [ [a, b, c, d, e, f, g, h, i] ] ; 

If [Test 56 I I D[Test, A] 9! 0, Coef f Ineq = Append [Coef f Ineq, Test > 0] ] ; 
, {a, 4}, {b, 4}, {c, 4}, {d, 4}, {e, 4}, {f, 4}, {g, 4}, {h, 4}, {i, 4}]; 
Reduce [Coefflneq] 

3 

Out[101]= A > - 
8 
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